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I .  INTRODUCTION 

The  singularity  expansion  method  (SEM)  has  received 
considerable  attention  by  workers  in  the  EMP  community.  This 
report  does  not  address  the  applications  of  the  method  to  EMP 
system  effects,  rather  we  elaborate  on  the  contribution  of  SEM 
to  an  increased  understanding  of  basic  scattering  theory. 

One  of  the  major  objectives  of  this  report  is  to  demon¬ 
strate  how  SEM  is  connected  to  general  scattering  concepts 
developed  by  workers  in  other  communities.  This  objective 
evolved  as  a  result  of  our  having  demonstrated  that  important 
theoretical  SEM  results  developed  for  electromagnetic 
scattering  were  obtainable  for  scalar  (acoustic)  scattering. 
Our  purpose  in  establishing  scalar  SEM  was  to  have  the 
benefit  of  dealing  with  scalar  scattering  to  address  some  of 
the  open  questions  that  exist  for  electromagnetic  SEM.  Some 
of  the  scalar  results  that  we  obtained  do  not  appear  in  the 
scalar  related  literature  we  utilized,  and  may  be  new. 

In  addition  to  our  scalar  SEM  results  as  well  as  the 
connection  to  related  theories,  this  report  presents  results 
from  two  separate  and  distinct  SEM  studies,  both  dealing 
specifically  with  electromagnetic  SEM.  One  study  was  an 
examination  of  the  relationship  between  class  1  and  class  2 
coupling  coefficients,  and  the  other  study  was  a  demonstration 
of  the  sensitivity  of  SEM  pole  location  values  as  a  function 
of  the  numerical  procedure  used  to  obtain  them. 

In  a  previous  work  (Ref.  1)  we  showed  that  class  2 
coupling  coefficients  yielded  results  that  were  in  conflict 
with  the  known  sphere  solution.  Despite  this,  it  seemed 
class  2  coupling  coefficients  might  still  be  of  some  use  if 
they  were  viewed  as  wrong  in  principle  but  possibly  a  good 
approximation  in  practice.  With  this  as  the  motivation  for 
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further  study,  we  obtained  the  following  results.  First,  we 
showed  that  after  a  time  which  is  the  sum  of  the  transit  time 
across  an  obstacle  of  an  incident  plane  wave  pulse,  plus  the 
time  corresponding  to  its  pulse  width,  class  1  and  class  2 
coupling  coefficients  yield  the  same  result.  Second,  for  the 
case  where  the  incident  field  was  a  step  function,  and  conse¬ 
quently  had  an  infinite  pulse  width,  we  showed  that  class  1 
coupling  coefficients  yield  the  correct  late  time  magneto¬ 
static  solution  while  class  2  coupling  coefficients  do  not. 

These  coupling  coefficient  efforts  are  presented  in  Appendix  D 
and  Appendix  E. 

We  now  return  to  our  scalar  SEM  investigation  and  describe 
it  in  more  detail.  It  is  our  view  that  scalar  SEM  theory 
would  be  comparably  established  relative  to  electromagnetic 
SEM  if  the  scalar  SEM  equivalent  of  the  following  electro¬ 
magnetic  SEM  results  could  be  obtained:  (1)  a  demonstration 
that  the  basic  concept  works  for  a  special  case;  the  work  of 
Baum  (Ref.  2),  which  initiated  electromagnetic  SEM  accom¬ 
plished  this  for  the  perfectly  conducting  sphere  scattering 
problem;  (2)  the  establishment  of  the  meromorphic  nature  of 
the  scattering  solution;  (3)  the  establishment  of  the  relation¬ 
ship  between  SEM  and  eigenmode  expansion  method  (EEM)  solutions 
in  order  to  obtain  explicit  relations  for  coupling  coefficients; 

(4)  a  demonstration  that  the  EEM-derived  solution  works  for  a 
particular  case;  and  (5)  the  establishment  that  cavity  SEM 
pole  locations  do  not  cause  nonphysical  results  for  the  external 
scattering  problem.  Results  (2) ,  (3) ,  and  (5)  appear  in  the  work 
of  Marin  and  Latham  (Ref.  3).  Later  Baum  elaborated  on  the 
relationship  between  EEM  and  SEM  in  Reference  4  and  obtained 
expressions  that  had  a  wider  use.  In  Reference  1,  we  established 
result  (4)  and  in  addition  we  further  clarified  the  relationship 
of  cavity  SEM  pole  locations  with  external  scattering  results. 

In  that  work  we  also  presented  a  more  useful  form  of  EEM  solutions 
than  had  previously  existed.  Much  of  the  remaining  literature 


concerning  electromagnetic  SEM  deals  with  numerical  calcul¬ 
ations  for  special  cases. 

In  this  report  we  emphasize  the  scalar  scattering  problem 
in  which  Neumann  boundary  conditions  are  imposed  on  the  surface 
of  the  scatterer.  We  choose  to  deal  with  this  scattering 
problem  because  it  directly  leads  to  a  scalar  integral  equation 
that  is  the  analogue  of  the  equation  that  provided  the  basis 
for  the  theoretical  investigations  contained  in  Reference  3, 
i.e.,  the  magnetic  field  integral  equation  (MFIE) .  As  will  be 
explained  shortly,  even  though  we  emphasize  the  Neumann  problem, 
the  structure  of  scalar  SEM  will  force  us  also  to  treat  the 
scalar  Dirichlet  scattering  problem  in  some  detail. 

Utilizing  both  the  Neumann  and  Dirichlet  scalar  integral 
equations,  we  establish  the  scalar  equivalent  of  results  (1)  , 
(3),  (4),  and  (5)  within  the  body  of  the  text.  Concerning 
result  (2) ,  the  meromorphic  character  of  the  scalar  scattering 
solution,  the  literature  dealing  with  scalar  scattering  treats 
this  feature  of  the  solution  as  a  foregone  conclusion.  A 
procedure  to  obtain  this  result  is  to  use  a  theorem  histori¬ 
cally  described  in  the  very  informative  review  article  by 
Dolph  and  Scott  (Ref.  5)  and  ultimately  attributed  to  Steinberg 
It  is  well  known  that  the  scalar  integral  equation  we  derive 
conforms  to  the  Steinberg  conditions  which  in  turn  guarantees 
the  desired  meromorphic  properties. 

Before  leaving  this  topic,  it  is  appropriate  to  describe 
some  of  the  work  in  more  detail.  To  provide  a  demonstration 
that  scalar  SEM  exists  for  a  particular  problem,  we  considered 
a  plane  wave  incident  on  a  spheroid  which  then  included  the 
sphere  as  a  special  case.  For  the  reasons  previously  discussed 
we  consider  the  case  where  Neumann  boundary  conditions  are 
satisfied  on  the  surface  of  the  spheroid.  We  show  that  the 
solution  obtained  by  the  EEM  is  the  same  as  the  standard 
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solution  one  obtains  by  separation  of  variables.  We  then 
specialize  this  solution  for  the  case  where  the  spheroid 
becomes  a  sphere.  We  rewrite  this  scalar  sphere  solution  in 
a  manner  which  exhibits  all  of  the  SEM  properties  that  Baum 
showed  for  the  electromagnetic  sphere  problem  in  Reference  2. 
Having  done  this,  we  were  immediately  in  a  position  to 
increase  our  knowledge  as  a  result  of  treating  the  scalar 
problem.  The  only  analytic  solution  for  scattering  from  a 
finite  object  that  it  is  possible  to  examine  in  the  electro¬ 
magnetic  case  is  the  solution  for  the  sphere.  For  the  sphere, 
we  find  that  the  eigenmodes  do  not  depend  on  frequency.  The 
scalar  sphere  EEM  solution  also  has  this  property;  however, 
the  scalar  spheroid  eigenmodes  do  depend  on  frequency. 

The  spheroid  EEM  solution  provided  information  in  another 
related  area.  Conditions  are  given  by  Ramm  in  his  very 
informative  review  article  (Ref.  6)  that  can  readily  be  inter¬ 
preted  to  be  a  set  of  sufficient  conditions  for  the  EEM 
solution  to  our  scalar  integral  equation  to  be  meaningful. 

Ramm  presents  enough  detail  in  that  article  for  us  to  conclude 
that  we  would  not  meet  the  described  sufficient  conditions 
unless  our  scalar  integral  operator  is  normal.  We  are  able  to 
show  that  this  is  the  case  when  the  scatterer  is  the  sphere, 
but  we  were  unable  to  show  this  when  the  object  was  the 
spheroid.  For  the  spheroid  we  show  that  the  integral  operator 
is  complex  symmetric  and,  despite  this,  the  EEM  solution  is 
shown  to  be  the  standard  separation  of  variables  solution. 

An  aspect  of  the  scalar  EEM  solution  worth  noting  is  the 
scalar  pseudosymmetric  analogue  of  the  electromagnetic 
pseudosymmetric  theory  presented  in  Reference  1.  We  show  that 
the  set  of  eigenvalues  for  the  Neumann  integral  operator  for 
the  external  scattering  problem  is  the  same  set  of  eigenvalues 
one  obtains  for  the  Dirichlet  integral  operator  for  the 
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interior  scattering  problem.  This  set  equality  implies  the 
existence  of  purely  imaginary  zeros  of  the  exterior  Neumann 
eigenvalues.  Just  as  in  electromagnetic  SEM,  these  extraneous 
zeros  must  be  shown  not  to  contribute  to  the  resulting  scalar 
SEM  expansion  and  we  provide  this  demonstration. 

We  now  present  a  more  detailed  description  of  how  we 
relate  scalar  SEM  to  general  scattering  concepts  that 
workers  in  other  communities  have  developed.  In  particular, 
we  will  connect  SEM  to  a  theory  of  scattering  developed  by 
Lax  and  Phillips  (Refs.  7  and  8).  This  theory  deals  with 
solutions  valid  over  a  volume  as  opposed  to  solutions  on  a 
surface  and  it  still  has  the  satisfying  quality  that  it  leads 
one  to  view  the  internal  scattering  (cavity)  problem  and  the 
external  scattering  problem  in  a  manner  that  exhibits 
considerable  uniformity.  A  connection  to  the  Lax-Phillips 
theory  acts  as  a  connection  to  a  chain  of  theories  as  effort 
on  their  part  was  made  to  connect  to  other  prominent 
scattering  theories  and  effort  on  the  part  of  workers  outside 
the  EMP  community  was  made  to  connect  to  the  Lax-Phillips 
theory. 

Our  introduction  to  the  Lax-Phillips  theory  came  from 
reading  the  review  article  by  Dolph  and  Scott  (Ref.  5) .  This 
informative  article  not  only  summarized  important  features 
of  the  Lax-Phillips  theory  but  connected  that  theory  to  other 
theories;  however,  it  only  minimally  mentioned  SEM.  We 
were  motivated  to  pursue  the  connection  because  we  had  been 
speculating  on  the  insights  into  SEM  one  could  obtain  if  a 
theory  such  as  that  of  Lax-Phillips  existed.  Lax  and  Phillips 
as  well  as  Dolph  and  Scott,  made  special  effort  to  relate  to 
an  approach  that  has  been  much  studied  in  quantum  mechanical 
scattering.  This  approach  utilizes  the  concept  of  a  scattering 
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operator  and  the  S-raatrix.  Another  work  that  makes  a  connec¬ 
tion  between  the  S-matrix  approach  and  the  Lax-Phillips  theory 
is  the  book  by  Nussenzveig  (Ref.  9).  This  book  contains 
detailed  material  concerning  applications  of  the  S-matrix 
approach  and  presents  a  general  review  of  other  efforts  where 
it  has  been  utilized. 

This  report  demonstrates  that  scalar  SEM  theory,  based  on 
the  EEM  expansion,  is  related  to  the  Lax-Phillips  theory.  In 
particular,  we  prove  that  the  set  of  complex  eigenvalues  that 
play  a  central  role  in  the  Lax-Phillips  theory  is  exactly  the 
same  set  as  the  one  consisting  of  the  (nonextraneous)  zeros 
of  the  eigenvalues  of  the  surface  integral  equation  with  the 
latter  set  being  the  SEM  pole  locations.  At  this  point,  one 
of  the  noteworthy  connections  between  Lax-Phillips  theory  and 
S-matrix  theory  is  that  the  complex  eigenvalues  have  been  shown 
to  be  in  one-to-one  correspondence  with  the  poles  of  the 
S-matrix.  We  have  now  established  that  SEM  pole  locations  are 
in  one-to-one  correspondence  with  these  S-matrix  poles. 

Another  important  feature  of  the  Lax-Phillips  theory  is  an 
asymptotic  expansion  (valid  for  late  times)  that  has  as  its 
only  time  dependence,  exponential  functions  having  an  argument 
that  is  the  product  of  the  complex  eigenvalues  and  the  time 
variable.  The  set  equality  we  proved  concerning  the  complex 
eigenvalues  and  SEM  pole  locations,  as  well  as  the  form  of  the 
SEM  expansion,  allow  us  to  conclude  a  term-by-term  equality 
of  the  two  expansions.  This  conclusion  is  based  on  the 
argument  that  the  common  time  functions  employed  in  both 
expansions  are  linearly  independent  functions. 

The  connection  between  the  Lax-Phillips  theory  and  scalar 
SEM  based  on  the  EEM  approach  benefits  both  efforts.  The 
Lax-Phillips  approach  has  a  more  advanced  theoretical  founda¬ 
tion.  As  an  example  of  this,  some  conditions  have  been 
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established  on  the  shape  of  scattering  surfaces  to  which  this 
theory  can  be  applied.  This  information  had  not  yet  been 
determined  by  the  SEM/EEM  approach.  In  addition,  quantative 
and  qualitative  estimating  techniques  have  been  developed  for 
the  complex  eigenvalues.  Finally,  as  has  already  been 
discussed,  the  Lax- Phillips  theory  has  been  connected  to  other 
scattering  theories.  The  SEM/EEM  approach  contributes  to  the 
scattering  problem  by  providing  explicit  expressions  for  the 
expansion  coefficients  and  these  expansion  coefficients,  as 
well  as  SEM  pole  locations,  have  been  numerically  determined 
by  workers  in  the  EMP  community.  More  generally,  workers 
in  the  EMP  community  have  developed  the  capability  to  numer¬ 
ically  obtain  SEM  solutions  for  scattering  shapes  that  are 
beyond  analytic  treatment. 

A  generalization  of  the  arguments  presented  for  the  scalar 
case  may  readily  be  possible  for  the  electromagnetic  case 
since  the  Lax-Phillips  theory  has  been  extended  to  electro¬ 
magnetic  scattering  (Ref.  10) .  We  would  be  inclined  to 
continue  focusing  attention  on  scalar  SEM  because  of  the 
untapped  information  contained  in  the  explicit  analytic  scalar 
solutions,  the  analogue  of  which  does  not  exist  for  electro¬ 
magnetic  scattering.  We  certainly  have  not  exhausted  the 
information  in  the  spheroid  solution  presented  in  this  report. 

Finally,  we  describe  a  result  that  we  were  able  to  obtain 
due  to  the  fact  that  we  have  made  the  described  connection 
between  scalar  SEM  and  Lax-Phillips  theory.  We  were  able  to 
obtain  a  formal  proof  that  the  SEM  poles  are  simple,  and  this 
has  long  been  identified  as  an  open  question  by  workers  in 
the  EMP  community.  This  formal  proof  is  totally  constructive 
in  nature  and  such  a  proof  can  alternately  be  thought  of  as 
providing  a  sufficient  set  of  conditions  for  the  desired 
conclusion.  The  sufficient  condition  nature  of  the  proof  is 
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implicit  in  that  certain  formal  manipulations  are  assumed 
valid  and  that  the  sufficient  conditions  are  those  implied  to 
make  the  manipulations  valid.  For  the  case  of  the  sphere, 
we  have  verified  that  these  manipulations  are  valid. 
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II.  A  SCALAR  THEORY  OP  THE  SINGULARITY 
EXPANSION  METHOD 

The  presentation  of  the  material  in  this  section  is  facil¬ 
itated  by  referring  to  Figure  1. 


Figure  1.  Separation  of  space  into  an  interior 
and  exterior  region. 

In  this  figure  we  introduce  a  surface,  S,  that  separates 
all  of  space  into  an  interior  region,  Vj,  and  an  exterior 
region,  VE.  At  this  point  S  is  just  a  mathematically  con¬ 
structed  surface;  however,  as  this  presentation  proceeds,  S 
will  correspond  to  a  physical  surface  on  which  boundary  con¬ 
ditions  are  satisfied  and  it  will  also  have  shape  requirements 

placed  on  it.  Also  in  Figure  1  are  sources  denoted  f„  and  fT, 

f  f  E  T 

which  are  nonzero  on  finite  volumes  VE  and  V^.  contained  within 
V£  and  Vj.  We  are  interested  in  finding  solutions  to  the 
scalar  wave  equation  in  each  region,  which  is  given  by 

“  A  - 2^-5 -  *  f  (r,t)  (1) 

a  c  at2  a  ~ 
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where 


a  *  E,I  for  reVE  or  reVj 

Employing  the  Laplace  transform  with  the  convention 


f  (s) 


f (t)e'st 


dt 


(2a) 


f  (t) 


(2b) 


where  Br  indicates  an  appropriate  Bromwich  path,  we  obtain 

(V2-y2)  ia(r,Y>  =  fa(£,Y)  (3) 

where  y  =  s/c.  In  the  remaining  portion  of  this  section  we 
will  omit  the  ~  notation.  We  now  introduce  the  Green's 
function  satisfying  the  equation 

<V2-y2)  G(r,r',Y)  -  -6  (r  -  r')  (4) 

which  is  given  by 


G(r,r* ,y) 


e“Y| r  -  r ' | 
4rr  |  r  -  r '  | 


(5) 


This  Green's  function  is  used  in  conjunction  with  the  identity 

7  .  [GV<f>  -*  VG]  »  GV2$  -  <P„V2G  (6)  \ 


Substituting  Equations  (3)  and  (4)  into  Equation  (6)  and 
integrating  over  Va  we  have 


[GV<|>  -  $  VG]  dV 


4>a  (r ' )  +  Qa(r') 


(7) 


where 

Qa  (r ’ )  =  J  fa(r)  G(r,r'  ,y)  dV  (8) 

vf 

a 

At  this  point  we  make  a  distinction  between  the  cases  a  =  E,I.' 
In  particular,  we  wish  to  apply  the  divergence  theorem  to  the 
left-hand  side  of  Equation  (7) .  In  order  to  do  this,  it  is 
necessary  to  identify  the  surfaces  bounding  the  volume  over 
which  the  integration  occurs.  For  the  integration  over  V^., 
the  bounding  surface  is  S  and  we  have 


[GV4>t  -  ^VGldV 


[GV$I  -  t^VGldS 


(9) 


For  the  integration  over  V£,  we  identify  the  bounding  surface 
as  being  S  as  well  as  a  spherical  surface  that  has  a  radius 
approaching  infinity.  We  now  have 


[GV<(>E  -  4>EVG]dV  = 


[GV4>e  -  <{>E VG ]  dS  +  IL 


(10) 


where 


I 


L 


ilim 

£-+CO 


[G V<|>E  -  4>EVG]  r2 


sine  d6d$ 


(11) 
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with  0,  <}),  and  r  being  spherical  coordinates  defined  in  a 
coordinate  system  having  its  origin  at  some  finite  distance 
from  S.  At  this  point  we  will  consider  that  the  surface  S 
is  the  physical  surface  corresponding  to  a  hard  acoustic 
scatterer  so  that  pj  and  <j>E  satisfy  the  Neumann  boundary 
conditions  on  S.  That  is. 


3  p 

“  *  v*a  =  IK  =  0  (12) 

In  addition,  <{>_  (r)  ,  in  the  terminology  of  Reference  7,  must 
E 

satisfy  the  "y-outgoing"  condition,  which  is  expressed  as 

<J>E  (r)  ~  K(e,4>)r-1e-Yr  (13) 

for  large  r.  For  a  y-outgoing  solution  for  <f>E  one  can  show 
that 


ar  •  [GV4>e  -  <4>eVG]  ~  F(r,r',y)e  2yr 
for  large  r  and  further  that  F  is  such  that 

lim  r2F(r,r',y)e'2Yr  =  0 

r  4oo 


(14) 


(15) 


for 


Rey  ^  o 


(16) 
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The  significance  of  Equation  (16)  is  that  it  defines  the 
region  for  a  Bromwich  path  that  in  turn  allows  us  to  conclude 
that 


JL  "  ° 


(17) 


It  is  worth  noting  that  the  specification  of  a  region  for  the 
Bromwich  path  has  always  been  inherent  in  the  more  usual 
Fourier  treatment  of  the  wave  equation.  The  Sommerfeld  radi¬ 
ation  condition  is  satisfied  if  Equation  (16)  is  satisfied,  as 
can  be  seen  by  making  the  substitution 


y  =  -ik 


(18) 


This  substitution  is  consistent  with  the  Fourier  transform 
convention 


f  (w) 


-l 


f (t) elwt  dt 


(19a) 


f  (t)  =  Jr  /  f  (u))e~lu,t  dco  (19b) 

yBr 

Both  Equations  (14)  and  (15)  are  valid  under  this  substitution 
and,  accordingly,  the  region  for  the  Bromwich  path  which  allows 
IL  to  vanish  is 

Jtmk  >  0  (20) 

For  y  or  k  restricted  to  the  identified  regions  so  that  IL  =  0, 
utilizing  the  Neumann  boundary  conditions  expressed  in 
Equation  (12) ,  and  changing  the  notation  so  that  the  roles  of 
r  and  r'  are  interchanged,  we  combine  Equations  (7) ,  (9) , 

(10) ,  and  (17)  to  obtain 


*0(r)  + 


/  ["c(E 


’)  •  V,G(r,r')]4)a(r,)dS'  =  -Qa(r)  (21) 


Equations  (8)  and  (21)  are  sufficient  to  give  (r) 
meaning  without  any  explicit  evaluation  of  integrals.  From 
Equation  (8)  we  see  that  Qa  (r)  is  independent  of  the  presence 
of  the  surface  S,  so  by  letting  the  surface  integral  in 
Equation  (20)  vanish,  we  conclude  that 


-Qa(r)  = 


.me,  \ 
a  <£> 


(22) 


inc 

where  <j>  (r)  is  the  field  radiated  by  the  sources  without  the 

surface  being  present.  Substituting  Equation  (22)  into 
Equation  (21)  yields 


L 


[na(r‘)  •  P'G(rv,r')]  ^(r')dS'  =  ^nc  (£v>  (23) 


and  we  have  changed  the  notation,  temporarily  replacing  r  by 
£v,  in  order  to  emphasize  that  Equation  (23)  was  derived 
with  r  .  ranging  over  the  volume  V  .  We  will  reserve  the 
notion  r  to  represent  an  observation  point  on  the  surface  S. 
Next  we  consider  the  limit  as  rv  approaches  the  surface,  i.e. 
r,  in  both  the  left-  and  right-hand  sides  of  Equation  (23) . 


Both  4>a( rv)  and  4>a 


me 


(r  )  are  continuous  functions  as  this 


limit  is  taken  for  values  of  rv  in  the  appropriate  so  that 
we  need  only  consider 


=  £im 


L 


tna <£'> 


V'G  (£v»£'  >  1 


<p  (r')dS' 


(24) 


18 


The  limit  described  in  Equation  (24)  is  a  frequently  treated 
limit  and  our  major  focus  in  this  report  is  to  present  a 
simple  bookkeeping  scheme  to  identify  critical  plus  and  minus 
signs  without  addressing  the  subtleties  of  this  limit.  Refer¬ 
ring  to  Figure  1  we  see  that  the  sign  of  the  dot  product  of 
the  vectors  rv~r*  and  n^Cr')  is  the  same  when  r  ,  whether 

a  =  E  or  a  =  I.  It  is  the  sign  of  this  dot  product  that 
determines  the  sign  multiplying  the  term  i  <J>a(r)  which  is  part 
of  the  expression  that  will  be  used  to  describe  the  limit. 
Since  the  sign  of  the  dot  product  is  independent  of  the  choice 
for  a,  the  sign  multiplying  the  i  (£)  term  in  the  resulting 
expression  is  also  independent  of  the  choice  for  a .  Using 
this  information  we  have 


I 


a 


2 


[  n  ( r ' ) 
a  —  * 


V' G (r ,  r '  )  ]<J>a(r'  )dS* 


(25) 


and  we  see  that  the  common  sign  is  a  negative  sign.  We  also 
note  that  the  factor  j  is  appropriate  only  if  r  does  not 
correspond  to  a  pathological  point  on  the  surface,  e.g.,  a 
tip  or  an  edge.  We  now  introduce  the  notation 

n(r')  =  nx  (r ' )  =  -nE(r')  (26) 

=  n(r  ’ )  •  V'G(r,r ' )  (27) 


K  $(r) 


ij>  (r '  )dS  ’ 


Combining  Equations  (23)  through  (28) ,  we  have 


(28) 
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(29) 


where 


le  4>e 
s  Y 


-  *E 
me 


L1  O1 


=  <$>* 

me 


(30) 


(31) 


(32) 


In  Equations  (29)  through  (32)  we  have  introduced  the  subscript 
"s"  to  the  operator  Ls  in  order  to  indicate  that  we  are  treat¬ 
ing  the  scalar  scattering  problem.  The  reason  for  this  addi¬ 
tional  notation  is  that  we  are  now  in  a  position  to  make  an 
important  distinction  between  the  scalar  scattering  problem 
and  the  vector  scattering  problem. 

From  Reference  1,  we  see  that  the  vector  perfect  conductor 
scattering  problem  treated  by  the  MFIE  has  a  generic  represen¬ 
tation  that  is  equivalent  to  Equations  (29)  and  (30).  It  is 
given  by 


and 


me 


(33) 


Tmc 

±1 


(34) 


(35) 


The  point  we  wish  to  make  is  independent  of  the  specific 
representation  of  l“  and  we  refer  the  reader  to  Reference  1 
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if  he  is  interested.  In  Equations  (33)  and  (34)  the 
quantities  J_,  JT,  jinc,  and  J^nc  have  analogous  meanings 
inc  inc 

to  4>e,  4>e  ,  and  4>  .  If  we  consider  that  the  same 

surface  depicted  in  Figure  1  was  perfectly  conducting  and 
fE  and  fj  were  replaced  by  electromagnetic  rigid  sources, 
fE  and  ,  then  we  can  describe  the  electromagnetic  quan¬ 
tities  as  follows.  The  quantity  JE  is  the  total  current 

density  induced  on  the  exterior  of  the  surface  that  had  as 

inc 

its  prime  excitation  fE,  and  JE  is  a  current  density  that 

would  be  excited  by  f£  at  the  same  mathematical  points 

describing  the  scattering  surface;  however,  the  physical 

surface  is  considered  to  be  removed.  The  quantity  is 

the  total  current  density  induced  on  the  interior  of  the 
inc 

surface  by  ,  which  again  is  the  current  density  excited 

by  fj.  on  the  mathematical  surface  S,  but  the  physical  surface 
is  considered  to  be  removed.  We  will  shortly  present  essen¬ 
tial  details  of  an  eigenmode  solution  to  Equation  (31)  as  we 
present  eigenmode  solutions  to  the  scalar  scattering  problem. 

We  now  consider  the  scalar  eigenvalue  equation 


La  <s>a  =  xa,s  <f>a 
s  n  n  n 


a  =  E,I 


(36) 


and  in  Appendix  A  we  define  an  inner  product  between  two 
functions  f  and  g  and  introduce  the  standard  notation  to  rep¬ 
resent  this  inner  product,  i.e.,  (f,g).  This  inner  product 

E  t 

is  used  in  the  determination  of  the  adjoint  operator  L 

S 

which  is  also  presented  in  Appendix  A.  The  eigenmode  expan¬ 
sion  method  (EEM)  solution  to  Equation  (36)  is  expressed  in 
terms  of  the  eigenmodes  and  eigenvalues  satisfying 

Equation  (36) ,  as  well  as  the  eigenmodes  satisfying  the 
adjoint  equation 
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Lat$at 
s  rn 


.  a,  s*  ,  at 

* 


a  =  E ,  I 


(37) 


The  corresponding  EEM^ solutions  are 


co 

V-  £ 

n=l 


\Tn  ,Tinc' 


E ,  s  n 
ln 


(38) 


<},1  = 


CO 

E 

n=l 


\Yn  me/ 


I,s 


I 

) 

n 


n 


(39) 


Using  results  from  Reference  1 ,  we  see  that  the  EEM 
solution  to  Equation  (31)  is 


(40) 


where  the  quantity  ijE,J.  >  is  a  pseudo-inner  product  which 

|  n  inc  j 

is  defined  in  Reference  1  and  its  explicit  representation  is 

not  required  for  us  to  make  our  desired  point.  The  eigen- 
F  E  V 

mode  and  eigenvalue  X^'  satisfy  the  equation 


E  S  =  E,v  E 
v— n  n  -n 


(41) 


where  LE  is  defined  in  Reference  1.  A  feature  that  makes 
the  connection  between  scalar  and  vector  EEM  as  well  as  SEM 
more  dramatic,  is  the  scalar  analogue  of  the  vector  pseudo- 
symmetric  theory,  presented  in  Reference  1.  In  order  to  make 
the  analogue,  we  retrieve  essential  results  from  Reference  1. 
First  we  introduce  the  vector  interior  equation  given  by 
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(42) 


l/j1  =  A  1 '  VJI 
V— n  n  — n 


To  make  the  desired  point  we  present  the  equation  satisfied  by 
~E 

the  Jn  which  appeared  in  the  pseudo-inner  product  contained  in 
the  EEM  solution  given  by  Equation  (40)  .  The  equation  satis¬ 
fied  by  this  quantity  was  obtained  from  adjointness  considera¬ 
tions  that  are  presented  in  detail  in  Reference  1 .  Retrieving 
the  desired  equation  from  that  reference,  we  have 


lJje 

V— n 


X 


E'viE 

n  — n 


(43) 


Comparing  Equations  (42)  and  (43),  we  conclude  the  set  equality 


|xE'Vl 

1  n  f 

l  n  ) 

(44) 


For  the  scalar  case,  there  is  also  a  connection  between 
the  interior  and  exterior  eigenvalues  that  appears  to  be  fun¬ 
damentally  different.  However,  viewing  the  SEM  consequences 
of  the  relationship,  shows  that  it  is,  in  fact,  fundamentally 
the  same.  In  Appendix  B,  we  present  the  details  to  support 
the  following  statements.  First,  we  cite  the  surface  eigen¬ 
value  equation  for  the  interior  Dirichlet  problem.  It  is 


,e+V  . 

s  n  n  n 


(45) 


where  is  defined  in  Appendix  B,  the  superscript  D  is 

affixed  to  the  eigenvalue  to  indicate  Dirichlet  boundary 

E  + 

conditions,  and  the  crucial  point  is  that  the  Ls  turns  out 
to  be  the  adjoint  of  the  Neumann  exterior  integral  operator. 
Comparing  the  complex  conjugate  of  this  equation  with 
Eguation  (37) ,  we  conclude  the  set  equality 
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ixI.s,D 
l  n 


»5'‘} 


(46) 


which  states  that  the  set  of  exterior  Neumann  surface  eigen¬ 
values  is  the  same  set  as  the  set  of  interior  Dirichlet 
surface  eigenvalues.  In  addition,  we  conclude  that  the  set 
of  exterior  adjoint  Neumann  eigenfunctions  are  the  complex 
conjugate  of  the  set  of  interior  Dirichlet  eigenfunctions. 
That  is , 


i  jEt*  ( 

)  I  i 

i^n  f 

(47) 

Using  this  result  we  can  write  Equation  (38)  in  a  manner 
analogous  to  Equation  (40)  as  follows: 


where 


E 

<? 


E 

■n  —  1 


J .1  ,E  I 

lJn,<?inci  x  E 


,E  } 
inc  ) 


R*>L)  * 


(48) 


(49) 


The  SEM  consequences  of  the  analogous  relations, 
Equations  (40)  and  (48),  will  now  be  discussed.  First  we 
express  the  following  relations  in  a  form  that  exhibits  a 
desired  amount  of  uniformity.  Equations  (38),  (39),  (40), 
and  (48)  can  be  written  as 


y«(s)  = 


oo 

E 

n=l 


N“(s) 

-^7T^n(s] 

A  n 


8  =  s ,  a  =  E, I 
8  =  v,  a  -  E 


(50) 
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n°(s)  =  Kt<=).*ino<s>)  «i»> 

Kv'sl  -  {S(s)'ii„cls)}  l51b) 

Nsls>  *  K(s)'*Lc<s>}  (51c> 

ysn(s)  =  4>n(s)  (51d) 

y“(s)  =  (j>a  (s)  ( 51e) 

s 

yLis>  m  Zn(s)  <51£) 

y®(s)  =  £E  (s)  (51g) 


The  argument  s  is  the  Laplace  transform  variable  while  the 
subscript  s  indicates  scalar.  We  arrive  at  the  desired 
SEM  related  residue  series  by  taking  the  inverse  of  that 
equation  using  the  Laplace  transform  convention  presented  in 
Equation  (2) .  From  independent  analyses  as  discussed  in  the 
Introduction,  we  know  that  y„(s)  is  a  meromorphic  function  of 

p 

s  and  from  additional  arguments  we  can  confine  the  poles  of 

y„(s)  to  the  left-half  s  plane  with  the  presence  or  absence 
p 

of  poles  lying  on  the  imaginary  axis  requiring  special  atten¬ 
tion.  The  argument  that  there  are  no  right-half  plane  poles 
is  related  to  causality.  It  is  argued  that  the  Bromwich  path 
can  be  closed  by  increasingly  large  right-half  plane  semi- 

Ot  3  £ 

circles  and  further  that  the  integral  of  yg(s)e  over  these 
semicircles  vanishes  for  all  t  less  than  the  time  when  the 
scattering  object  is  initially  struck  by  an  incident  field. 
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This  allows  us  to  close  the  Bromwich  path  to  the  right  for 
these  values  of  time.  Because  of  causality,  we  do  not  wish  to 
enclose  any  poles  since  they  v.  .uld  yield  residue  contributions 
to  the  inverse  transform  for  nonphysically  interpretable 
values  of  time.  There  is  a  connection  between  causality  and 
the  radiation  condition  in  that  the  radiation  condition 
arguments  led  to  the  specification  that  the  Bromwich  path  lie 
to  the  right  of  the  singularities,  thus  making  the  zero 
residue  response  possible. 

We  now  employ  analytic  continuation  and  close  the  Bromwich 
path  to  the  left  for  a  range  of  positive  times.  The  results 
for  this  case  are  formal  in  that  we  will  simply  make  any 
mathematical  assumption  required  to  obtain  the  results  desired. 
If  we  proceed  along  this  line,  a  question  arises  as  to  the 
utility  of  our  results.  In  a  subsequent  section  we  will  make 
some  of  these  results  rigorous,  and  in  addition  we  will  explain 
how  these  results  are  consistent  with  those  obtained  by  anal¬ 
yses  which  did  not  make  these  assumptions.  These  results  will 
also  be  shown  to  be  in  agreement  with  a  known  sphere  solution. 

To  obtain  the  desired  form,  we  look  for  the  totality  of  the 

ot  8 

zeros  of  the  eigenvalues  X  ,R(s)  and  we  denote  an  arbitrary 

a,  B,T  .  n 

zero  as  s',  where 


(s 


a ,  6  ,  T 
nn' 


) 


0 


n'  =1,2, 


.  NT(n) 


(52) 


and  the  superscript  T  is  used  to  indicate  totality.  As 
a  result  of  Equation  (44)  for  the  vector  case  and  Equation  (46; 
for  the  scalar  case,  we  expect  a  decomposition  of  the  total 
set  of  zeros  as  follows 
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{  sa'M 

U 

(sa,e,EXj 

l  nn'  ) 

l  nn'  ) 

l  nn’  f 

(53) 


The  zeros  s^,  are  the  significant  zeros  and  the  zeros  having 
the  superscript  EX  attached  are  the  extraneous  zeros.  For 
the  exterior  problem,  the  significant  zeros  are  the  ones  having 
a  negative  imaginary  part  and  the  extraneous  ones  are  purely 
imaginary.  For  the  interior  problem,  the  significant  ones  are 
purely  imaginary  and  the  extraneous  ones  are  in  the  left-half 
plane.  Equations  (44)  and  (46)  demonstrate  the  existence  of 
the  extraneous  zeros  for  the  exterior  problem  because  we  know 
that  the  interior  surface  eigenvalue  contained  in  each  of 
these  equations  has  purely  imaginary  zeros.  The  value  of  the 
pseudosymmetric  argument  presented  in  Reference  1  and  in 
Appendix  B  is  to  eliminate  nonphysical  residue  results  due 
to  the  extraneous  zeros  for  the  exterior  problem.  More 
specifically,  in  Reference  1  it  is  shown  that 


E,v,EX 

■nn' 


0 


and  in  Appendix  B,  we  show  that 


E ,  s  ,  EX. 
s  ,  ) 

nn 


0 


(54) 


(55) 


Assuming  that  the  significant  zero  is  simple,  we  write  the 
residue  arising  in  the  inverse  Laplace  transform  corresponding 


.a,  S 


a,  8  a 


to  the  significant  zero  s  ,  as  ,(r)e 

5  nn'  nn,J 6nn' 


sa,8t 
snn ' 


where 


_a,  6 
ann  * 


N 


a 

6 


(s 


a,  6v 
nn' ; 


X' 


a  >  6 

n 


(s 


a,  6 

nn ' 


) 


(56) 
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(57) 


a  i  \  a  ,  a ,  B  . 

y6nn'(£>  =  ySn  snn 1 


.a,  6  f  cl  r  3 » 
Xn  snn 1 


dA  a '  8 ( s ) 
n 

ds 


.  _  _  —CL  t  B 

is  -  snn 1 


(58) 


Using  these  residues,  we  write  the  response  as 


y“(r,t)  =  s“(r,t)  +  R“(r,t) 


(59) 


where 


s“(r,t) 


N(n)  sa,Bt 

a*'?y“  ,  (r)  e  nn * 

nn’-^enn’  — 


EE 

n=l  n'=l 


(60) 


and  Ro(r,t)  is  just  the  remaining  function  needed  to  augment 

^  Ct  OL 

the  sum  Sg(r,t)  in  order  to  represent  y„(r,t).  The  real 
assumption  is  that  Sg(r,t)  converges  in  some  meaningful  sense 
and  then  it  also  dominates  Rg(r,t).  By  comparing  the  surface 
approach  used  in  the  EMP  community  to  the  volume  approach 
used  by  workers  outisde  of  that  community  (Ref.  8) ,  we  can 
infer  that  Sg(r,t)  does  have  the  desired  property.  Finally, 
we  note  that  an  alternate  form  for  the  scalar  expansion,  one 
that  bears  a  closer  resemblance  to  the  more  prominent  electro¬ 
magnetic  SEM  expansion,  is  presented  in  Appendix  C. 
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III.  BACKGROUND  MATERIAL  FOR  THE  CONNECTION 
TO  THE  LAX- PHILLIPS  THEORY 


To  facilitate  the  desired  connection  between  the  theory 
just  presented  and  the  Lax-Phillips  theory  which  is  a  volume 
approach,  we  introduce  the  volume  eigenvalue  equations  for 
both  the  interior  and  exterior  scalar  scattering  problems 


( V2  -  y  2  )  <j>  =0 

an  an 


a  =  E,  I 


(61) 


For  the  interior  problem,  a  =  I,  the  boundary  condition 


3  d> 

an 

3n 


0 


(62) 


on  the  surface  as  well  as  certain  volume  behavior  require¬ 
ments,  e.g.,  requirements  which  force  us  to  reject  the  explicit 
solutions  obtainable  for  separable  coordinates  that  become 
unbounded,  are  known  to  lead  to  denumerable  sets  of  eigen¬ 
functions  UIn}  and  {yIn}-  Furthermore,  for  the  eigenfunction 
equation  convention  used  in  Equation  (61) ,  it  is  also  known 
that  the  YIn's  are  purely  imaginary  and  occur  in  complex  con¬ 
jugate  pairs.  To  clarify  this  and  to  provide  a  point  of 
reference  for  a  subsequent  discussion,  we  note  that  the  more 
common  form  of  the  interior  eigenvalue  equation  is 

,v2  +  kln>tn  '  0  1631 

which,  when  the  same  boundary  conditions  and  volume  behavior 
requirements  are  imposed,  it  is  known  to  lead  to  a  discrete 
set  of  eigenfunctions  and  eigenvalues  with  the  further  result 
that  the  k_  '  s  lie  on  the  real  axis. 
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Returning  to  Equation  (61),  we  consider  the  exterior 
scattering  problem.  The  fact  that,  subject  to  appropriate 
boundary  conditions,  there  exists  only  a  denumerable  set  of 
eigenfunctions  and  eigenvalues  is  not  as  well  known  to  be  the 
case  for  the  exterior  problem  as  it  is  for  the  interior 

problem.  The  boundary  conditions  that  lead  to  these  denu¬ 
de, 


merable  sets  are 


En 


8n 


=  0  on  the  surface  the  the  y-outgoing 


condition  expressed  as 


KEn 


-1  _YEnr 
Kn(0,<j,)r  e  bn 


(64) 


for  large  r.  As  discussed  in  Reference  8,  we  also  have  the 
additional  information  that 


ReyE  <  0  (65) 

'n 

The  requirement  on  the  shape  of  the  scatterer  is  an  issue 
that  has  received  attention.  Many  of  the  cited  properties  have 
been  proved  when  the  object  is  star  shaped,  i.e.,  a  point 
within  the  object  can  be  found  from  which  a  straight  line  can 
be  drawn  that  connects  this  point  to  any  other  point  within 
the  volume  bounded  by  the  surface  of  scatterer.  Star-shaped 
surfaces  include  convex  surfaces.  There  has  been  some  work  and 
some  conjecture  concerning  the  applicability  of  the  work  to 
surfaces  described  as  confining  or  nonconfining  (Ref.  8).  It 
should  be  noted  that  the  predominant  situation  of  an  imperfectly 
sealed  enclosure,  i.e.,  a  finitely  thick-walled  enclosure  con¬ 
taining  an  aperture,  is  not  a  star-shaped  surface. 

We  now  restrict  our  attention  to  surfaces  for  which  we 
have  the  desired  discrete  spectrum  for  the  exterior  scattering 
situation  and  utilize  the  identity 
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(66) 


V 


[<j>  VG  -  GV<j> 
an 


an 


<p  V2G-G  V2<J> 
an  an 


Substituting  Equations  (4)  and  (61)  into  Equation  (66)  ,  we  have 


V*  [<f>  VG-GV<J>  ]  =  (Y2  -  Y2  )4>  G  -  6  (r-r‘ )  <p 

an  an  an  an  - Tan 


(67) 


First  we  consider  this  equation  for  the  interior  problem  and 
integrate  both  sides  over  the  interior  volume  and  use  the 
divergence  theorem  as  well  as  the  boundary  condition  given 
by  Equation  (62)  to  obtain 


(r') 


VG(r,r')<j)In(r)dS 


<t>,  GdV 
In 

(68) 


Interchanging  the  notation  r  and  r'  and  taking  the  limit  as 
r  approaches  the  surface,  we  have 


In 


In 


GdV’ 


(69) 


where  L*  is  defined  by  Equations  (27),  (28),  (31),  and  (32). 

We  now  rewrite  Equations  (36)  and  (69)  emphasizing  the  Y 
dependence.  Equation  (36),  after  changing  the  dummy  index 
to  m,  is 


(70) 


and  Equation  (69)  is 


I.'<Y)*In<r)  = 


<1Xn'y2)II(Y'YIn-£) 


(71) 
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where 


II(’f’TIn'El 


<t>In(r  '  )G(y,  |  r-r '  |)dV' 


(72) 


Rather  than  discuss  the  connection  between  SEM  and  the  volume 
approach  at  this  stage,  i.e.,  having  only  provided  material 
for  the  connection  to  the  interior  problem,  we  will  obtain 
an  equation  that  has  the  same  significance  for  the  exterior 
scattering  problem  as  does  Equation  (71)  for  the  interior 
problem. 

To  obtain  the  desired  equation,  we  consider,  two  mathemat¬ 
ical  surfaces  to  be  introduced  in  the  exterior  domain.  One 
is  the  smallest  sphere  that  can  circumscribe  the  scattering 
object,  and  the  other  is  a  sphere  whose  radius  we  will  allow 
to  approach  infinity.  Initially,  we  will  not  make  use  of  the 
smaller  sphere  and  will  integrate  both  sides  of  Equation  (67) 
over  the  volume,  VE  ,  which  is  bounded  by  the  scattering 
object  and  the  large  sphere,  denoted  Sr  ,  and  having  a  radius 
r.  Performing  this  integration  we  obtain 


♦En^  + 


/Ss<E>  ' 


VG(r,r')<t.En(r)  dS 


L 


11,2  -  L  *EnGdV  -  Ir 

VE 


(73) 


where 


I 


r 


ilim 

£*-►00 


[GV*En“  VG5dS 


(74) 
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Now  the  "y-outgoing"  condition  leads  to  conclusions  analogous 
to  those  expressed  by  Equations  (14-16).  That  is,  for  large  r 


[Gv*En  -^VG] 


En 


Fn(r,r* ,y, 


YEn)e 


~(Y 


+  YEn)r 


(75) 


where  is  such  that 
n 


dim 

r->oo 


r  Fn(r,r*,Y,YEn)e 


-(Y+YEn)r 


=  0 


(76) 


for 


Recognizing  the  Rey., 

c*n 


Ee(y+YEn) 
=  -UeyEn|, 
Rey  2  | Rey 


>  0 

this  condition  becomes 
En  ^ 


(77a) 


(77b) 


Since  lReYEnl  is  a  nondecreasing  function  of  n,  we  see  that 
the  specification  of  the  Bromwich  path  moves  to  the  right  as 
n  increases.  Despite  this,  we  can  make  the  choice  of  the  path 
independent  of  n,  at  least  for  some  time,  by  employing  the 
following  argument.  We  can  analytically  continue  any  function 
of  y  that  is  of  interest  back  to  the  line  Rey  =  —  |  Re'YEi  1  where 
|ReYEil  is  the  smallest  I ReYEn I  for  aH  n*  can  perform 
this  analytic  continuation  in  the  finite  y  plane  without 
encountering  any  singularities  since  all  the  Y_  have  negative 
real  parts.  We  now  see  that  with  the  Bromwich  path  initially 
chosen  according  to  Equation  (77) ,  it  follows  that 

Ir  =  0  (78) 
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We  now  turn  our  attention  to  the  volume  integral  appear¬ 
ing  in  Equation  (73)  and  denote  this  integral  as  i”.  We 
write  this  integral  as  the  sum  of  three  terms  and  in  order 
to  do  this,  we  utilize  the  previously  defined  smallest 
circumscribing  sphere.  That  is,  we  express 


I 


n 

E 


<j>En(r)G(Y,  jr-r*  |  )  dV 


alternately  as  the  sum  of  the  terms 


(79) 


where 


rn 


(80) 


=  Him 

r~>-oo 


L  »EnGdV 
va 

f  <*E„G-VdV 

vb  . 


(81) 


(82) 


Him 

r~>-co 


(83) 


To  complete  the  definitions  expressed  by  Equations  (80) ,  (81) , 

and  (82) ,  it  is  necessary  to  define  V  ,  V.  ,  and  f  .  The 

a  D  ..  n 

quantity  V  is  the  volume  between  the  surface  of  the  scatterer 

a 

and  the  circumscribing  sphere  whose  radius  we  denote  as  a. 

The  quantity  is  the  volume  between  the  circumscribing  sphere 
and  the  increasingly  large  sphere  having  the  radius  r.  The 
quantity  f  is  a  function  that  is  constructed  as  the  product 


of  the  functions  describing  the  large  r  asymptotic  behavior 
of  <}>En  and  G.  For  large  r,  the  asymptotic  behavior  of  4>En  is 
given  in  Equation  (64)  and  we  just  note  that  the  Kn  given  in 
that  equation  is  a  function  of  angular  variables  which  will 
now  be  defined.  Let  us  consider  a  spherical  coordinate  system 
defined  as  having  its  origin  at  the  center  of  the  circum¬ 
scribing  sphere  and  let  us  denote  the  usual  angular  variables 
as  6  and  <j>.  Then  we  can  write 


Jlim 

r-*-°° 


sine 


f 


n 


(84) 


For  our  purposes  it  is  sufficient  to  represent  f  as 


f 


n 


*n(r',e,*,y,YEn> 


+  Y)r" 


(85) 


and  we  shall  see  shortly  that  the  exact  form  of  is  not 
important.  Substituting  Equation  (85)  into  Equation  (84) , 
we  have 


Lbf 


=  Ve' 


’YEn^ In 


sing 


(86) 


where 


A 


n 


sinG 


(87)  i 

i 

i 


and 


i 
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[sin9  =  Urn  f 

n  r+»  -4 


r  -(Y+YFn)r" 
e  En  dr" 


=  ilim 

jr-^oo 


~(Y+YEn)a  '(^En)r 
e  -  e 


(88) 


If  Equation  (77)  is  satisfied,  then 


sing  ,  ^  ,-l  -(Y+YEn)a 

*n  "  ,Y+YEn>  e 


(89) 


Equation  (73)  can  now  be  written  as 


^En (£*  }  +  /  nE(r)  *  VG(r,r')'J>En(r')dS 

s 


“  (Y2-YEn)Ic(Y'YEn^')  +  ^Y  ~  YEn)  XD  (Y '  YEn '  — ' }  (90) 


where 


=  I2(E’'Y'YEnl 


Y'W 


(91) 


and 


4  ■ 


An'E' 


'Y’YEn,e 


-<Y+YEn)a 


(92) 


The  significance  of  the  terms  contained  in  Equations  (91)  and 
(92)  is  that  they  are  not  singular  when  y  =  YEn  >  and  this  is 
still  the  case  when  r’  approaches  the  surface  of  the  scatterer. 
We  are  now  concerned  with  taking  this  limit  in  Equation  (90) 
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and  we  introduce  the  same  notation  change  that  was  used  in 
taking  the  same  limit  of  Equation  (68) .  That  is,  we  inter¬ 
change  r  and  r*  notation,  use  Equation  (26) ,  as  well  as  the 
limiting  arguments  preceding  that  equation  to  obtain 


Ls(Y)<t,En(^  =  +  (Y"YEn)IS{Y'YEn'^) 


where  L 


E 

s 


is  defined  in  Equation 


(31)  . 
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IV.  CONCLUDING  ARGUMENTS  CONCERNING  THE  CONNECTION 
BETWEEN  THE  TWO  THEORIES 


At  this  point  we  will  retrieve  appropriate  equations  from 
the  preceding  section  and  rewrite  them  in  a  more  convenient 
form.  We  begin  by  presenting  Equation  (70),  without  modifi¬ 
cation,  to  make  it  convenient  for  reference. 


L“(y)<t>“(Y,r)  =  x“'S(Y)4»°‘(Y,r)  a  =  E,I  (70) 

o  in  in  m  ‘ 

We  rewrite  Equations  (71)  and  (93)  as 

LS(Y)*an(r)  =  (Y-Yan)Fan(Y,r)  (94a) 

and 

Lg  (Y) =  {Y+YIn)Fjn( Yf£)  (94b) 

where 


FIn  =  -(Y+YIn)II(Y'YIn^) 


Fln  =  {Ym-Y)II(Y'YIn^) 


FEn  =  ( Y+YEn]  Fc  ( Y '  YEn }  +IS(Y'YEn'^ 


'•  A  major  conclusion  that  we  wish  to  draw  concerns  the 
following  sets.  One  set  is  described  following  Equation 

sa'® 

in  conjunction  with  the  definition  Y^f  =  nJ>—  ,  and  we 
write  this  set  as 


(95a) 


(95b) 

(96) 


(53) 


The  other  set  is 


{ y  \  =  (y  I  (V2-y  2)<J>  =  0  +  bo^arY  \  (98) 

I  am)  l  am  1  am  am  conditions) 

The  first  conclusion  that  we  wish  to  draw  is  the  set  equality 

(99) 

To  prove  that 

y  l  iC'i  <i°°> 

we  consider  the  interior  problem  first.  We  do  this  to  empha¬ 
size  a  major  difference  between  the  interior  and  exterior 
problem.  This  difference  manifested  itself  in  the  need  to 
present  Equation  (94b)  for  the  interior  problem.  Consider  a 
typical  YIm  and  let  y  =  YIm  in  Equation  (94a)  to  obtain 

Ls(Wim<£>  -  0  noi) 

Letting  y  =  in  Equation  (94b)  yields 


“  0  ll02> 

I  s 

Equation  (101)  implies  that  there  exists  some  A £  (y)  such 
that  =  0#  while  Equation  (102)  implies  that  there 

X  s  X  s 

exists  some  Aj'  (y)  such  that  Xj'  (~YIm)  =  0.  This  not  only 
proves  Equation  (100)  for  a  =  I,  but  it  also  proves  that  if 


then 


v  i 

*  Im 


/  I ,  s  I 


nn  ’ 


(103a) 


(103b) 


The  conclusion  expressed  by  Equation  (103)  is  acceptable  for 
a  =  I  but  would  not  be  acceptable  for  a  =  E.  The  reason  for 
making  this  observation  is  to  emphasize  the  need  for  the 
bounding  arguments  presented  in  the  previous  section  that  led 
to  Equation  (94a)  for  a  =  E,  and  at  the  same  time  eliminated 
the  need  for  an  equation  corresponding  to  Equation  (94b)  for 
the  exterior  problem.  Substituting  y  =  y  into  Equation 
(94a)  leads  to 


L^Yj^Jr)  =  0  (104) 

S  Em  Em  — 

which  (by  exactly  the  same  arguments  applied  to  Equation  (101)) 
leads  to  the  result 


w 


(105) 


however,  there  is  no  result  corresponding  to  Equation  (103b). 
This  is  necessary  in  that 


Rey„  <  0  (106) 

cj  m 

while,  as  is  well  known  for  the  interior  problem, 

ReyIm  =  0  (107) 


Equation  (107),  together  with  the  fact  expressed  in  Equations 
(103a)  and  (103b),  simply  states  that  y,  occur  in  complex 


conjugate  pairs.  If  an  equation  corresponding  to  Equation 

(103b)  were  true  for  the  exterior  problem,  then  because  of 

Equation  (106),  we  would  be  led  to  contradictions  concerning 

the  existence  of  right  half-plane  values  of  YEm*  The  bounding 

arguments  presented  in  the  previous  section  precluded  this 

undesirable  conclusion.  Finally,  we  note  that  the  more  familiar 

form  of  Equation  (107)  is  obtained  by  using  the  complex  plane 

rotation  described  by  Equation  (13) ,  which  was  presented  in 

conjunction  with  the  complex  Fourier  transform.  The  statement 

equivalent  to  Equation  (107),  after  making  the  indicated 

substitution  k_  =  iy _  is 
im  Im 


ImkT  =  0 
Im 


(108) 


Equation  (108)  is  the  more  familiar  result  corresponding  to 
internal  resonances. 

Equations  (103)  and  (105)  prove  Equation  (100) .  To  com¬ 
plete  the  proof  of  Equation  (99) ,  we  will  prove  that 


M  -kk 


(109) 


To  accomplish  this,  we  form  the  following  function 


.a,  x  i  a,s,  .at*.  a,s 

l>  (r  )  =  /  G(r  ,r' ,y  , )  (y  ,,r')dS' 
n  — v  J  — v  —  nn'  n  nn'  — 


(110) 


where  *  indicates  complex  conjugate  and  <t>  (y,r' )  satisfies 

the  adjoint  eigenvalue  equation.  The  adjoint  operator  is 


given  in  Appendix  A,  and  the  adjoint  eigenvalue  equation  is 
given  as 


T  at ,  .  , at  ,  .  .  a , 

Ls  (Y)*n  (Y'^  =  V 


(y)  (y  ,r) 


(111) 


As  in  the  previous  section,  subscript  v  is  attached  to  the 
r  to  indicate  that  r  varies  over  the  appropriate  three- 

Q, 

dimensional  volumes  VT  or  V„.  Direct  substitution  of  ip  (r  ) 

IE  n  —  v 

given  by  Equation  (110)  into 


2  a  u 

(V  - y  7  )  P  (r  ) 
v  nn  rn  — v 


(112) 


shows  that  this  equation  is  satisfied  for  r  in  the  appropriate 

£  V 

open  volumes  V£  or  .  In  addition,  <p  (£v)  can  readily  be 
seen  to  satisfy  the  Y-outgoing  condition.  For  Equation  (105) 
to  be  true,  it  remains  to  show  that  the  Neumann  boundary  condi¬ 
tions  on  S  are  satisfied.  To  show  this  we  take  the  gradient 
of  both  sides  of  Equation  (110)  to  obtain 


V  ip 
v 


01  (r  )  =  f 

n  ~v  Je 


/  C-(r  ,  r  1 
v  -v  — 


ra'f) 
nn  ’ 


.at* 


(Y 


a ,  s 


nn  — 


r  1  )  dS  1 


(113) 


Identifying  r'-as  a  point  on  S  which  will  be  approached  by  r  , 
and  defining  na(r)  as  the  appropriate  normal  defined  at  the 
point  r,  consistent  with  the  convention  depicted  in  Figure  1, 
we  have 


Aim  n  (r)  •  Vv<j>“(rv)  =  Aim 

r  -  r  a  r 

— v  —  — v  - 


n  fn(  r ) 

r  J  a 


V  G ( r  ,r 1 
v  —  v  — 


,y 


nn 


?) 


at* 

n 


(ya,s 
Ynn 1  ' 


r ’ )dS' 


(114) 
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The  appropriate  1/2  factors  come  from  the  limit  expressed  by 
the  right-hand  side  of  Equation  (114)  so  that  the  following 
equation  is  true: 


~  ,  a  ,  ,  _  at*  ,  ci.s,  .at*,  a ,  s  , 

n  ( r )  *  (r)  =  L  (y  . )  <p  (y  ,  ,r) 

a  —  n  —  s  'nn'  Tn  '  'nn'  — 


(115) 


First  we  recognize  that  the  left-hand  side  of  Equation  (115) 

is  the  desired  normal  derivative.  Next  we  substitute  y=ya,f 

nn' 

into  Equation  (111)  and  then  take  the  complex  conjugate  of 
the  resulting  equation  to  obtain 


Tat*,  a,s.  ,  at*  ,  a,s  .  .a,s,  a.s,  ,at* ,  a.s  . 

Ls  Ynn'  -An  ^nn'  •£>  (ll6) 


Combining  Equations  (115)  and  (116)  ,  we  obtain 


3iJja{r)  4.* 

n  —  ,ct,s.  a,s. ,at*  a,s  x 

-35—  -  X„  aOT,)*n  (Ynn..E> 


(117) 


where  we  have  used  the  convention  that 


3<l>“(r) 

3n 


-n  (r)  •  Vtp  (r) 
a  —  n  — 


(118) 


Referring  to  the  definition  of  the  y^f  given  in  Equation  (97)  , 
we  see  from  Equation  (117)  that 


3<J^(r) 

3n 


=  0 


(119) 


Equations  (112)  and  (119)  yield  the  desired  result 
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(120) 


which  in  turn  implies  that  Equation  (109)  is  true.  Equa¬ 
tions  (100)  and  (109)  together 'prove  the  desired  result 
expressed  by  Equation  (99) . 

Shortly,  we  will  interpret  the  significance  of  Equation 
(99) .  However,  we  will  first  point  out  the  intrinsic  con¬ 
sistency  of  the  equations  that  led  to  this  result.  There 
were  two  intermediate  results  that  occurred  during  the  course 
of  this  study  that  appeared  as  detriments  to  drawing  desired 
conclusions.  One  was  the  fact  that  the  integral  given 
by  Equation  (79)  became  unbounded.  The  other  initially 
appeared  as  a  nuisance  factor  when  we  determined  that  the 
operators  turned  out  not  to  be  self-adjoint  as  described 
in  Appendix  A.  The  analysis  presented  in  this  section  showed 
that  both  of  these  apparently  undesirable  results  turned  out 
to  be  very  desirable.  The  unboundedness  of  In  is  the  under¬ 
lying  reason  that  precluded  the  consideration  of  -y  as  be- 

g  g  Juill 

longing  to  W^n' '*  The  nonself -adjointness  of  l“  greatly 
facilitated  the  proof  of  Equation  (109) . 

Returning  to  the  interpretation  of  Equation  (99) ,  we  first 
prove  that  the  elements  of  this  single  set  which  we  now  denote 
as  {y“}  for  notational  convenience,  i.e.. 


(a  \ 

1  ) 

i  a ,  s\ 

<Y_  ( 

— 

1  Y  (  = 

<Y  ,  > 

(  am) 

nn' 

(X 

contain  the  pole  locations  of  the  <p  that  satisfy  Equations 
(29)  and  (30),  which  we  rewrite  as 


L“(YHa(Y,r)  =  4“nc(Y,r) 


(122) 
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Prom  arguments  already  discussed,  we  know  that  4>  (y,r)  is  a 
memomorphic  function  of  y  having  all  of  its  poles  in  the 
closed  left-half  y  plane.  First  we  will  consider  the  possi¬ 
bility  that  a  pole  at  the  location  pna  has  multiplicity  M  . 


Next  we  multiply  both  sides  of  Equation  (122)  by  (y-p  ) 


M 


.na 


and  then  set  y=  p 


na 


Assuming  we  have  not  chosen  an  incident 


field  having  the  same  pole  location,  the  resulting  equation  is 


La(p  )Ra(r)  =  0 
s  'na'  n  — 


(123) 


where 


M 

_a ,  .  ,  .  na  ,  a  ,  , 

Rn(r)  =  (Y~Pna)  4>  <Y,r) 


(124) 


Y  =  P 


na 


Equation  (123)  implies 


P  e 
'na. 


(125) 


which,  according  to  Equation  (121)  ,  implies 


P  e 
*■  na 


(126) 


At  this  stage  we  have  not  presented  an  argument  that  all  of 

ct 

the  members  of  (Yn)  are  pole  locations. 

A  major  goal  of  this  work  is  to  facilitate  the  use  of 
complementary  efforts  related  to  SEM  that  have  been  generated 
by  different  communities.  In  Reference  8,  the  following 
asymptotic  expansion,  expressed  in  the  notation  of  this  work, 
is  given  as 


(127) 


4>e  (r ,  -t )  ~  S£(r,t) 

where 


SE(r,t) 


CO 


m=l 


^  c r ) 

n,m  — 


(128) 


and  t  =  ct.  4>-p(r,t)  is  said  to  behave  asymptotically  like 

hi  — 

S£(r,t)  for  large  t. 

After  performing  only  a  cursory  search,  we  found  no 
explicit  representations  for  the  c  1 s  nor  did  we  find  esti¬ 
mates  on  errors  introduced  by  terminating  the  sum  as  a 
function  of  either  r  or  t.  Despite  our  inability  to  obtain 
this  information,  the  existence  of  Equation  (127)  together 
with  our  result,  Equation  (99) ,  benefits  our  state  of  knowl- 

E 

edge.  A  formal  comparison  of  S  (r,t),  given  by  Equation  (60), 
after  the  following  substitution, 


_  E ,  s^_  E ,  s 

S'  ,  t  y  ',T  Y„  T 

nn  nn  Em 

e  =  e  =  e 


(129) 


with  S£(r,t)  given  by  Equation  (128)  yields  the  following 
formal  conclusions.  The  sum  of  all  of  the  terms  multiplying 
YEmT 

a  particular  e  ,  i.e.,  allowing  for  degeneracy  of  the 
complex  eigenvalue  YEm»  in  the  representation  S£(r,t)  given 
by  Equation  (128),  is  equal  to  the  sum  of  the  terms  multi- 

y  x 

plying  the  same  e  m  in  the  representation  of  S^(r,t)  given 
by  Equation  (60).  This  formal  equality  has  mutual  benefits. 

K 

The  terms  contained  in  the  representation  S  (r,t)  have 
explicit  representations  in  terms  of  surface  quantities  and 
furthermore  they  are  numerically  calculable. 


46 


The  benefit  in  going  the  other  direction  is  that  we  now 
have  additional  reasons  to  believe  that  the  remainder  term 

g 

R  (r,t)  becomes  unimportant  for  large  t.  This  result  has 
s 

frequently  been  assumed  to  be  true  by  workers  in  the  EMP 

E 

community,  based  on  the  argument  that  Rv(r,t)  was  the  inverse 
transform  of  an  entire  function.  This  argument  suffers  in  that 
a  very  simple  entire  function  can  be  identified  whose  inverse 

Oi 

transform  would  dominate  Sg(r,t)  in  a  controlled  manner.  This 
function  for  the  scalar  case  is  simply 


fE  (Y)  =  Ae 


(130) 


If  f  (y)  were  the  additive  entire  function  for  the  case 

where  the  time  dependence  of  the  incident  field  had  a  Dirac 

delta  time  dependence,  then  the  inverse  transform  of  f  (y) 

E 

convolved  with  the  time  dependence  of  a  general  incident  field 

(X 

could  be  chosen  to  yield  R^(r,t)  which  is  larger  than  the 
s“(r,t).  All  that  would  be  required  would  be  to  choose  the  con¬ 
stants  A  and  T  so  that  this  occurred.  The  asymptotic  nature  of 
SE(r,t)  precludes  this  possibility  for  the  scalar  scattering 

problem  and  suggests  there  is  no  hidden  surprise  in  the 
E 

Rv(r,t)  for  the  vector  scattering  problem. 

Before  leaving  the  formal  connection  between  S_(r,t)  and 

E  " 

S  (r,t),  two  additional  comments  are  in  order.  The  first  is 

s  ~ 

just  to  note  that  the  SE(r,t),  which  results  from  the  volume 
approach,  contains  the  features  discussed  in  the  Introduction. 
The  second  is  that  the  form  of  S_  is  very  indicative  that 
<J>„(r,y)  has  only  simple  poles  at  the  location  y_  .  This 

£j  —  EIT1 

follows  from  taking  a  term-by-term  Laplace  transform  of  the 
sum  S_(r,t).  We  will  shortly  present  a  separate  structured 
argument  that  this  is  actually  the  case. 


Our  argument  that  4>e(£,y)  has  only  simple  poles  requires 
a  number  of  clearly  identified  mathematical  statements  to  be 
true.  A  crucial  mathematical  requirement  is  that  Equation 
(29)  have  the  eigenmode  solution  expressed  by  Equation  (38) . 
Changing  the  notation  of  Equation  (25) ,  we  ask  whether  a 
solution  ^(y,r)  exists  to  the 


Lf(Y)!HY,r)  =  F(y,r> 

b 


(131) 


for  a  given  F(y,r)  and  furthermore  we  ask  that  the  solution 
can  be  expressed  as 


4*  (Y/£) 


E 


(132) 


The  conditions  that  tj;  exist  for  a  specified  F  can  be  found  in 
discussions  of  Fredholm  operators  since  Lg(y)  is  such  an 
operator.  The  conditions  for  the  validity  of  the  EEM  expressed 
by  Equation  (132)  require  further  discussion.  A  recently  dis¬ 
tributed  report  by  Ramm  (Ref.  6)  presented  conditions  that  an 

F 

operator  very  closely  related  to  Lg(y)  could  be  inverted 
according  to  Equation  (132).  In  particular,  he  treated  the 
operator  presented  in  Appendix  A.  He  concluded  that,  for 
a  suitable  class  of  scattering  surface  shapes,  the  root 
system  of  K  formed  an  appropriate  expansion  basis.  He  then 
emphasized  the  limitations  of  this  conclusion  in  that  the  root 
system  consists  of  both  eigenvectors  and  root  vectors.  In 
order  to  utilize  a  standard  EEM  expansion,  it  is  necessary 
that  there  be  only  eigenvectors  and  no  root  vectors  in  the 

root  system.  His  conclusions  concerning  the  root  system  for 

•f*  FI 

K  can  be  extended  to  apply  to  the  operators  L  and  Li  and  we 

S  b 
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now  focus  on  L  .  In  order  that  Equation  (132)  be  valid,  it  is 
S  £ 

necessary  that  the  root  system  of  Lg  contain  only  eigenvectors 
and  no  root  vectors.  Ramm  (Ref.  6)  gives  a  set  of  sufficient 
conditions  for  this  to  be  the  case,  and  shows  that  all  of  these 
are  readily  satisfied  with  the  exception  of  the  condition  that 
Lg  is  normal.  Ramm  points  out  that  whether  this  sufficient  condi¬ 
tion  is  satisfied  depends  on  the  shape  of  the  scatterer.  He  cites 
some  of  his  earlier  work  which  demonstrated  for  the  analogous 
Dirichlet  problem  that  a  spherical  scatterer  yielded  a  normal 
operator.  Later  in  this  work,  we  examine  whether  Lg  is  normal 
for  both  the  sphere  and  spheroid.  For  the  sphere,  we  conclude 
that  it  is  normal  while  for  the  spheroid  we  could  only  establish 
that  it  is  a  complex  symmetric  operator.  Despite  this,  we 
conclude  that  the  EEM  solution  for  the  spheroid  yields  the 
correct  solution  that  is  obtained  by  the  standard  separation  of 
variables  procedure. 

We  now  consider  Equation  (132)  as  applied  to  a  special  case 
of  Equation  (131).  This  special  case  is  Equation  (94a)  with 
a  -  E  and  the  dummy  index  changed  to  i.  For  this  situation 
Equation  (132)  becomes 


♦ei'E*  “ 


V  Y~Y£i. 

»!!'•<« 


(*^+<T'E> 'FEi <Y.£>)  (133) 


We  have  already  shown  that  for  some  n,  say  N, 


A5'S(*Bi>  =  0 


(134) 


In  order  that  the  right-hand  side  remain  finite  and  nonzero 

E  s 

at  Y  =  ye^  ,  we  conclude  that  A  '  (y)  has  only  a  simple  zero 
at  yE.. 


Referring  to  Equation  (132)  for  a  source  F  (y,r)  cor¬ 
responding  to  an  incident  field  have  a  Dirac  delta  time 


a 


dependence,  we  make  the  following  observation  concerning  the 
meromorphic  function  ^D(Y,r).  The  poles  of  >i'D(Y,£)  corre¬ 
sponding  to  the  zeros  of  the  eigenvalues  are  simple  poles.  It 
is  conjectured  that  for  many  surfaces  these  are  the  only  poles 
For  the  case  where  the  incident  field  has  a  more  general  time 
dependence,  we  simply  convolve  the  inverse  transform  of  the 
Dirac  solution  4'D(Y»r)  with  the  specified  time  dependence. 

We  conclude  this  section  by  making  an  important  distinc¬ 
tion  between  the  vector  and  scalar  scattering  problems.  We 

have  just  presented  sufficient  conditions  for  the  eigenvalues 
E  s 

A  '  (Y)  to  have  simple  zeros.  Equation  (60),  which  represents 

a  standard  SEM  expansion  for  all  cases  of  interest,  vector  and 

scalar,  interior  and  exterior,  is  based  on  the  assumption  that 

all  of  the  eigenvalues  An'p(y)  have  simple  zeros.  In  this  sec 

tion  we  have  presented  sufficient  conditions  for  this  to  be 
E  s 

the  case  for  An'  (y)  .  Identically,  the  same  arguments  could 

I  s 

be  applied  to  conclude  that  An'  (y)  had  only  simple  poles. 

All  that  would  be  required  would  be  to  use  Equation  (94a)  and 
(94b)  for  a  *  I  to  obtain  an  equation  analogous  to  Equation 
(133)  for  the  interior  problem.  It  now  remains  to  discuss 
the  vector  problem.  As  mentioned  in  the  Introduction,  there 
has  been  work  directed  at  extending  the  scalar  volume  approach 
to  the  vector  (electromagnetic)  case.  We  have  not  pursued 
that  vector  extension  to  determine  whether  the  vector  analogue 
of  Equation  (61)  has  been  established  for  the  exterior  vector 
problem.  If  it  has,  the  arguments  presented  in  this  work 
would  require  minimum  modification  to  draw  analogous  conclu¬ 
sions  to  the  ones  drawn  for  the  scalar  case.  The  primary 
modifications  would  be  to  utilize  vector  Green’s  theorems 
and  the  free  space  dyadic  Green's  functions.  This  might  be 
readily  accomplished;  however,  it  is  our  view  that  further 
insight  into  SEM  can  be  obtained  by  continuing  the  study  of 
the  scalar  problem.  We  have  two  underlying  reasons  that  led 
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us  to  this  view.  One  is  that  only  sufficient  conditions 
have  been  identified  that  allow  desired  SEM  conclusions  to 
be  drawn.  The  other  is  that  there  are  more  finite  shapes  for 
which  analytic  solutions  exist  for  the  scalar  scattering 
problem,  the  spheroid  is  only  one  of  them.  In  this  report, 
we  have  gained  some  information  by  examining  the  scalar 
spheroid  scattering  solution;  however,  we  have  not  exhausted 
the  information  contained  in  that  solution. 


\ 
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V.  PROLATE  SPHEROID 


In  this  section  we  consider  an  acoustically  hard  prolate 
spheroid  and  determine  the  eigenfunctions  and  eigenvalues  of 
the  surface  integral  operator.  These  eigenfunctions  have  an 
explicit  dependence  on  the  Fourier  complex  variable  k  (or  y 
if  Laplace  transforms  are  employed)  wnereas,  as  in  the  vector 
case,  the  eigenfunctions  for  the  special  case  of  a  sphere  are 
independent  of  k. 

The  scalar  integral  equation  for  an  acoustically  hard  body 
is 

<J>inC  (r )  ,  reS  (135) 


\  4>(r) 


/< 


3  . 


an' 


G ( r ' ; r )  dS'  = 


where 


G 


ikR 
e _ 

4tiR 


,  R  = 


I  £“£ 


9G 
3n ' 


n(r’)  •  V ' G 


and  n(r' )  is  the  outward  unit  normal  at  r ' .  The  corresponding 
eigenvalue  problem  has  the  form 

L4>  =  A  (136) 

n  n  n 


where 

L»n  =  (|  -  K)  »n 

K,n  =/  5FdS  <1371 

s 
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The  prolate  spheroidal  curvilinear  coordinates  5,  n ,  $  are 
related  to  the  Cartesian  coordinates  x,  y,  z  by  the  transformation 


I  1/2 


cos 


y  -  |  [(«2-l)  (l-n2)_ 


1/2 


sin  (J) 


z  =  2  dCn 


where  £  is  the  radial  coordinate  (1<^<'»)  ,  n(-l£n_^l)  and 
4)  ( 0_< 4>_< 2tt )  are  the  angular  coordinates  and  d  is  the  interfocal 
distance.  The  surface  of  the  spheroid  is  defined  by  {;  =  £.. 
The  correspondence  with  the  spherical  coordinates  r,  0  is 


£d  -*■  r,  n  -*■  cos  0  as  z,  -*• 


i.e.,  we  let  the  interfocal  distance  go  to  zero  and  keep  £d 
finite.  We  also  present  the  formulas  for  the  metric  coeffi¬ 
cients  hr,  h  , 

£  n 


h,  because  we  will  use  them  shortly: 
<t> 


h„  = 


•*( 


,2  2 

^  ~'i 

[2-l 


1/2 


h  = 
n 


/  2  2\1/2 
■  s  (») 


h.  = 


[{l-n2)  (r.2-l)] 


1/2 
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The  scalar  Hermholtz  equation 

(v2+k2)  $  =  0 


has  the  following  solution  for  outgoing  waves 


<5 


S 


mn 


(c,  n) 


r{3) 

mn 


(c,  £) 


( cosmij)  | 
|  sinmij)  | 


where 

c  =  j  kd 


(3) 

and  S  (c.n),  R  (c,£)  are  the  angular  and  radial  functions 
mn  mn 

respectively.  The  equations  they  satisfy  are  given  in  the 

monograph  by  C.  Flammer  (Ref.  11) .  In  the  limit  of  a  sphere, 

S  (c,n)  reduces  to  the  associated  Legendre  Polynominal  Pm(cos0) 
mn  , o \  / 1  \  n 

and  R  (£)  to  the  spherical  Bessel  Function  h  (r) .  The 
mn  n 

angular  functions  smn(c»ri)  satisfy  the  following  orthogonal¬ 
ity  relationship 


(c,  n) 


sm„ , (c, n) 
mn 


dn 


nn 


N 

mn 


(138) 


where  N  is  the  normalization  factor  given  in  Reference  11. 
mn 

In  terms  of  the  spheroidal  functions  the  Green's  function  G 
has  the  following  expansion 

eikR  ik  (2-6  ) 

TiF  =2^  spe(c'n)  spii(c'n,)  c°sp(H') 

p=0  £=p 


(c,C  ) 


(c,£) 


£  >  V  (139) 
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where  (c,£)  corresponds  to  jJ?(kr),  (c,£)  to  h^  (kr) 

and  N  ^  is  the  same  normalization  factor  as  in  Equation  (138) 

We  are  now  in  a  position  to  show  that  the  eigenfunctions 
of  L  defined  by  Equation  (137)  are: 


<J 


mn 


=  S  (c  n)  \cosm(P' 
mn  '  n  |  sinm'j) 1 


(140) 


Equation  (140)  shows  the  explicit  dependence  of  $mn  on  k  or 
Y(c=  (1/2)  kd  =  (i/2)  yd).  We  begin  by  noting  that  the  inte¬ 
gral  operator  involves  the  normal  derivative  of  the  Green's 
function  G 
given  by  Equation 


In  order  to  be  able  to  employ  the  representation 


;±39)  and  operate  on  (c,£r)  inside  the 


sum,  we  recall  that 


1 

2 


<&__ (r) 
mn  — 


9G 
9n  1 


dS’ 


S 


=  <5> 

mn 


(r) 


9G 
3n ' 


dS' 


(141) 


where  G(C1»n*  #4>*  in  the  second  line  has  £  +  i.e., 

its  unprimed  radius  vector  corresponds  to  a  point  off  the 
surface.  In  the  first  line  the  limiting  process  has  been 
completed  and  both  r  and  r1  correspond  to  points  on  the  sur¬ 
face  £  =  £^.  In  view  of  Equation  (139)  we  can  write 


55 


3G  _  1  8G 

3nT  ,  3T 


5'  =  5- 


CO  oo 


ik(2"6op) 

>  f~t  ~^v  sp^(c,n)  V(c'n,)  cosp^-^') 


p=0  J6=p 


he  dCx  pH 


Ml) 


Rpi'  «="«) 


€  >  5, 


If  we  now  notice  that 


dS  1  =  h  dn  1  h  ,  ,  d<}>  ’ 

n  4> 


h  .  h ,  ,  _ 

n1  4>  _  f2  . 


and  take  into  account  the  orthogonality  relationships  of  £ 
and  of  the  trigonometric  functions,  we  find  that 


lim 

5-5 


J 


S 

ikd  ir.2 


$  <r’)  dS* 

mn  —  dn 


I5!-1)  R™’Ml)  d?7  Rmi’(=-«l)  ♦„„<£> 


where  $  is  given  by  Equation  (140).  Combining  Equation: 
mn 

(137),  (141),  and  (143)  we  obtain 


mn 


- 1  -  ic(«H  CMJ  a^’Mi) 


(142) 


mn 

(143) 
;  (136), 

(144) 
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where  the  relationship  c  =  (1/2) kd  has  been  used.  In  order 

to  simplify  Equation  (144)  we  employ  the  equations  satisfied 

by  the  radial  functions  R^ ^ \  (Ref.  11)  alonq  with  their 

mn  mn  - 

asymptotic  expressions  for  large  £  (also  given  in  Reference  11) 
to  obtain  the  Wronskian  relationship 


dR(1)  ...  dR(3) 

(3)  mn  _  R(l)  mn  _  _  l 

mn  d£  mn  d£  ~  c  U2.]} 


(145) 


Combining  Equations  (144)  and  (145)  we  obtain 


X 


mn 


R(1) 

mn 


d£. 


R 


(3) 


mn 


Mi) 


(146) 


If  the  Laplace  transform  variable  y  is  used  we  simply  replace 
c  by  iyd/2  wherever  c  appears. 

Now  that  we  have  determined  the  eigenfunctions  and  eigen¬ 
values  for  Equation  (136)  by  utilizing  Equation  (142) ,  we  can  go 
back  and  give  an  explicit  expansion  of  9G/9n'  when  both  r 
and  r'  are  on  the  surface  of  the  spheroid.  The  answer  is 


3G 
9n ' 


®  1  _ 

ZEE 

£=0  p=0  o 


1-2  X 


2N 


El 


$> 


r) 


p£ 


p£a  v- 


$  (r '  ) 
p£o  -  ' 


(147) 


where  o  stands  for  odd,  even.  in  Appendix  A  we  show  that  for 

both  the  spheroid  and  the  sphere  the  adjoint  operator  L+  is 

★ 

equal  to  L  where  the  asterisk  signifies  complex  conjugation. 
For  the  sphere  we  also  show  that  L  is  normal.  In  general,  in 
order  to  solve  Equation  (135)  we  need  the  adjoint  eigenfunc¬ 
tions.  A  formula  for  the  adjoint  operator  is  derived  in 
Appendix  A. 
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In  the  limit  of  a  sphere  we  easily  obtain 


$ 

mn 


=  Pm(cos6) 
n 


i  cosm<}>  | 
I  sinm4>  j 


(148) 


?  ( i  \  1 

X  =  -  i(ka)  j  (ka)  h'1'  (ka) 
n  n  n 

=  -  (ya)2  in (ya)  k^(ya)  (149) 

where  f' (x)  i  df/dx.  The  eigenfunctions  given  by  Equations  (140) 
and  (148)  form  complete  orthogonal  sets  in  their  respective 
domains  (-l_<n<_l ,  0<$<2ir  for  the  surface  of  the  prolate  spheroid 
and  Oj<0<7r,  0<b<2n  for  the  sphere)  and,  consequently,  they  repre¬ 
sent  the  complete  solutions  to  the  eigenvalue  problems.  Along 
with  the  eigenvalues  they  can  be  used  to  solve  the  integral 


i 

where 

($r$inc)  =  f  *£(r)  $inc  (r)  dndcf) 

s 

and  $A,  X A ,  are  the  eigenfunctions,  eigenvalues  and  normal¬ 
ization  factors,  respectively.  Employing  Equations  (140)  and 
(146)  ,  we  can  show  that  the  scalar  field  evaluated  at  the 

surface  of  a  hard  spheroid  excited  by  an  incident  plane  wave 
ikz 

e  is  given  by  the  well-known  formula 


(yMC) , 
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c<q  -  X)  n^O 


£  Non  td^l  Ron)(c,Cl)]  Son(c'1)  Son(c'n) 


(151) 


by  noting  that  on  the  surface  of  the  spheroid 


’  L  iCT  son<c'1'  RoJ’ '='«!» 


and  that 


eikZsondn  -  2i%„«c'i»  Roi,|c'Ei> 


For  the  sphere  it  is  easy  to  derive  the  SEM  expansion  by 
first  employing  Equations  (148)  and  (149)  and  assuming  an 
incident  plane  wave  e  YZ 


(  P™(cos6)  ,e  Yz  c<?sm<}>\ 

*(£,Y)  --  -2 - S-n  ‘  ( cos  8 )  c?s  m<)> . 

iWa-(Ya)2in(Y«)k;(Ya)M 


Following  the  methods  employed  in  References  1  and  2,  we  arrive 
at  the  following  SEM  expansion 


4>(r,y)  =  e 


-yet 


(2n+l) (-1) 

F  1  /  y  ,  a 

(  nn1 


)  C'-vf  » 


P  (COS0) 


(152) 
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where 


Fix)  i  eVk'(x),  F'(x)  ,  g 

Ynn'  are  the  roots  of  kn(Ya)  =  0  and  tQ(=  -a/c)  is  the  instant 
at  which  the  wavefront  first  hits  the  sphere.  Equation  (152)  » 

corresponds  to  class  1  coupling  coefficients  as  expected  from 
the  vector  case . 
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VI.  NUMERICAL  RESULTS 


In  this  section  we  present  numerical  results  for  the 

zero  locations  corresponding  to  the  vector  sphere  problem. 

These  zeros  were  found  by  employing  the  method  described 

in  Reference  1.  Figure  2  shows  the  division  of  the  surface 

of  the  sphere  into  zones,  where  NL0NG  is  the  number  of  slices 

in  the  azimuthal  direction  and  NTTim  is  the  number  of  slices 

LAT 

in  the  polar  direction. 

Table  1  presents  the  numerical  results.  Notice  that  the 

search  routine  (Ref.  1)  locates  clusters  of  zeros  corresponding 

to  a  true  single  zero.  For  example,  for  n  =  1  and  NT =  4, 

LONG 

NT„m  =  2,  we  obtain  the  three  zeros  -1 . 0202+i0 . 0862 ,  -1.1012 
LAT  ' 

-i0.0522,  -0.9408-i0.0365,  whereas  the  real  zero  is  s^g  =  -1+iO. 
The  reason  for  the  clustering  is  that  zoning  transforms  the  sur¬ 
face  of  a  sphere  into  a  different  surface  and  the  three-fold 
degeneracy  (n=l,  m=-l,0,l)  is  resolved.  As  NL0NG  “*  °°»  NLAT  “ 
all  three  zeros  will  coalesce  into  the  one  zero  s11Q  if  we 
ignore  numerical  roundoff  error.  Also  notice  that  as  we  increase 
the  number  of  zones  a)  more  zeros  belonging  to  a  cluster  are 
found  and  b)  the  norm  of  the  error  of  the  average  cluster  value 
for  a  zero  decreases. 

The  rectangular  area  we  employed  is  defined  by: 

Coordinates  of  the  lower  right  corner:  (-1.55,  -0.30) 
Coordinates  of  the  upper  right  corner:  (0.00,  1.30) 
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I 


Figure  2. 


Zoning 

N 

LONG 


scheme.  For  the  figure 
8(M=45°),  =  2  ( A9=90° 


). 


1 

‘4 


TABLE  1.  POLE  LOCATIONS  AS  A  FUNCTION  OF  ZONING 


Designation 
of  zeros 


Real 

part  *(-l) 

..  . 

Imaginary 

part 

0.8787 

0 

0.5554 

0.8763 

0.9763 

0 

0.9306 

0 

0.5701 

0.8663 

0.4525 

0.8250 

1.3252 

0 

1.0013 

0 

0.9620 

0 

1.4398 

0.8410 

1.1146 

1.0637 

0.5374 

0.8897 

0.5128 

0.8523 

1.0202 

0.0862 

1.1012 

-0.0522 

0.9408 

-0.0365 

1.4621 

0.8491 

1.4610 

0.8483 

0.56912 

0.912328 

0.56913 

0.912322 

0.4892 

0.8586 

1.5479 

0.0033 

0.9813 

0.0023 

1.0396 

0.0103 

1.0269 

-0.0125 

1.4957 

0.8339 

1.4588 

0.8742 

0.5670 

0.8737 

TABLE  1 .  CONCLUDED 


Number  of  zones 

Zeros 

No.  of 

Designation 

Real 

Imaginary  i 

nlong 

nlat 

zeros 

of  zeros 

part  *(-l) 

part  | 

S21 1 

0.5016 

0.8648 

■■ 

s211 

,  ■ 

0.9242 

8 

iiv, 

8 

sno 

0.01258 

sno 

0.9921 

-0.0059 

sHHfpi 

sno 

1.0179 

-0.0067 

1 

s  1 21 

1.5063 

0.9291 

1 

S1 21 

1.4548 

0.8388 

C 

211 

0.5600 

0.9043 

s211 

0.5071 

0.8196 

■ 

'211 

0.4651 

0.8979 

s110  =  1,  s211  =  -0.5000  +  iO . 8660 ,  s121  =  -1.5000  +  i0.8660 
S220  =  _1-6961  +  iO.0000 
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APPENDIX  A.  ADJOINTNESS  RELATIONSHIPS 


In  this  appendix  we  derive  an  expression  for  the  adjoint 
to  the  scalar  integral  operator  corresponding  to  an  acousti¬ 
cally  hard  body  and  we  specialize  it  for  the  case  of  the 
spheroid  and  sphere. 

We  start  with  the  exterior  equation 


1 

2 


$  (r) 


/ 


3G  (r r) 

$  (r"  )  - - -  dS1 

3n 1 


S 


where  I  is  the  identity  operator. 

The  adjoint  is  defined  through  the  relationship 


(Al) 


(A2) 


where  we  have  used  the  inner  product  definition 


( f  t  g) 


g(r)dS 


i.e.,  we  can  define 


To  find  we  write 


( 5  K  $ ) 


K$<r)dS 


(r)  $  {r') 


S  S 


3G(r’;r) 

-  dS'dS 

3n ' 


(A3) 
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3G(r;r') 

ip*  (r ’  )  *(r)  — 5 -  dSdS  1 

—  —  dn 


/r  3G*(r.r')  V 
*(r' )  - -  dS' 

S 


where  the  penultimate  step  involved  an  interchange  between 
r  and  r ' . 

Thus,  according  to  Equation  (A2) 


=  J  \p  (r ' 


3G*(r  r’  ) 

>  — 3^r — ds' 


where 


3G (r .  r 1 )  3G(r'.r) 


i -  =  n  (r)  •  VG 


-  -(s +  Y) 


4iR~  ^*n(^ 


3G (r' .r) 


=  n(r' ) •  V’G 


=  (s  +  *) 


WT  5- ■><£’) 


The  interior  integral  equation  is 


i-  $(r)  +f  <t  (r ' 


3G (r ' ; r) 


dS'  =  (|  +  K)  $  E  L] 


i"  + 

and  thus  the  adjoint  operator  is  (1/2)  +  K  where  K  is  given 
by  Equation  (A4) . 


For  a  spheroid  one  can  use  Equation  (147)  to  cast  the  oper¬ 
ator  in  the  form 

L *-!*-/’  «<*•>£  T5T1  V£l  %<£’>  d— 1 

J  ~r  m 

m 

and  employ  the  procedure  outlined  in  this  appendix  to  show 
that 


t  * 

The  adjoint  eigenfunctions  of  L  are  then  For  a  sphere, 

in  particular,  one  can  choose  real  eigenfunctions  <j>m  and 

consequently  the  adjoint  eigenfunctions  are  identical  to 

the  4>  *s;  i.e.,  the  operator  L  is  normal, 
m 
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APPENDIX  B.  SCALAR  ANALOGUE  OF  THE  ELECTROMAGNETIC 
PSEUDOSYMMETRIC  THEORY 


We  begin  this  appendix  by  viewing  certain  properties  of 
the  eigenvalue  for  the  exterior  spheroid  scattering  problem, 
which  was  derived  in  Section  V.  This  eigenvalue  is  repro¬ 
duced  here  for  convenient  reference. 


W*’  -  <=-«!»  afr*™ 


(Bl) 


where 


c  =  |  yd  (B2) 

and  the  associated  functions  and  parameters  have  been  defined 
in  Section  V.  Our  effort  in  this  appendix  is  motivated  by 
examining  the  two  ways  in  which  the  zeros  of  ^mn(Y)  are 
obtained.  They  are 


d 


R<3> 

mn 


(c ,  = 


0 


(B3) 


and 


C  «o«l»  -  0  (B4> 

An  examination  of  whether  A  had  an  as  yet  unexplained  zero 

mn 

at  y  =  0  was  conducted.  We  found  that  the  multiplicative  y 
that  appears  in  this  eigenvalue  is  cancelled  by  an  aporo- 
priate  term  that  is  contained  within  the  dR^n  (c,C^)/d^^ 
factor.  Equation  (B3)  is  an  explicit  example  of  the  con¬ 
nection  we  proved  between  the  significant  zeros  of  the 
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eigenvalue  and  the  Lax-Phillips  complex  eigenvalues.  Equa¬ 
tion  (B4)  is  the  clue  to  pursue  the  results  developed  in  this 
appendix.  That  is,  it  yields  the  interior  volume  eigenvalues 
associated  with  the  Dirichlet  problem. 

We  now  begin  our  investigation  of  this  issue  by  deriving 
the  interior  Dirichlet  integral  equation.  First  we  combine 
Equations  (7) ,  (9) ,  and  (22) ,  employ  the  change  of  variables 
previously  described  in  which  r'  becomes  the  integration 
variable  and  r  an  observation  variable  in  V  ,  and  employ  the 

v  d  1 

Dirichlet  condition  ^  =  0  on  S,  to  obtain 

4>°<rv)  +  J G(rv,r*)oI(r')dS'  =  4>Jnc  frv>  (B5) 

S 

where  we  have  affixed  the  superscript  D  to  indicate  Dirichlet 
boundary  conditions  are  employed  and  we  introduce  the 
definition 


aI(r')  =  -njSrMV'^tr')  (B6) 

We  now  take  the  gradient  of  both  sides  of  Equation  (B5)  to 
obtain 

Vv<(>°(rv)  +  y"vv G(rv,r,)oI(r')dS'  =  Vv<^nc(rv)  (B7) 
S 

Identifying  a  point  r  on  the  surface  as  the  point  approached 
by  rv  ,  and  defining  n^.  (r)  at  the  point  r  ,  we  first  take 
-rij  (r)  dotted  into  the  preceding  equation  and  then  take  the 
limit  as  rv  approaches  r  to  obtain 
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1 

2 


o1  (r) 


3G(r,r' ) 
3n 


o1 (r ' ) 


=  Fa<^ 


(B8) 


where  we  have  used  Equation  (26)  and  the  definitions 


and 


n(r)  •  VvG(rv,r' ) 


SG( r,r') 
3n 


-“<£)  ’  Vi"°(£*>  ’  Fo(£> 


(B9) 


(BIO) 


The  sign  of  the  j  a  factor  in  the  described  limit  is  deter¬ 
mined  in  the  same  manner  that  was  described  in  treating 
Equation  (24) .  We  will  now  write  Equation  (B9)  as 


ho”1  * 


(Bll) 


and  this  equation  serves  as  the  definition  of  L  .  We  now 

Et 

compare  L^D  with  L  given  in  Appendix  A  and  find  that 


L 


ID 


(B12) 


We  now  consider  the  interior  Dirichlet  eigenvalue  equation 


LID°n 


I  _  XI,S,D„I 


n 


n 


(B13) 


Next,  we  take  the  complex  conjugate  of  Equation  (B13) ,  substi¬ 
tute  Equation  (B12)  into  the  result  to  find 
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(B14) 


L 


E+  I* 
s  an 


=  A 


I,s,D*  I* 
n  n 


Comparing  Equation  (B14)  and  Equation  (37)  ,  we  conclude  that 


and 


Ia1'3'0! 

\  n  / 

1  n  / 

(B15) 


{€*}  -  K}  (Bi6> 


Equation  (B15)  is  consistent  with  the  spheroid  eigenvalue 
result  expressed  in  Equation  (B4)  and  in  addition,  it  is  very 
significant  by  itself,  as  discussed  in  Section  II.  In  par¬ 
ticular,  it  implies  the  existence  of  the  extraneous  zeros  of 
the  external  eigenvalue,  which  we  will  now  argue  do  not 
contribute  to  the  desired  SEM  expansions. 

Specifically,  for  the  external  Neumann  surface  problem, 
we  want  to  show  that 


E(  E,s,EX) 
s  nn ' 


0 


(B17) 


£  g  EX  E  s 

where  s'.'  is  an  extraneous  zero  of  A  '  corresponding  to 
nn  n 

the  zero  of  the  interior  Dirichlet  problem.  To  accomplish 
this,  we  first  employ  the  identity 


V 


v  <t> 


E 

inc 


] 


0?  V2*1D  -  <pIDV2<t>E 

Tmc  n  n  me 


fB18) 


E  ID 

This  identity  is  true  for  general  functions  <f> .  and  <}>  ,  but 

p  me  n 

we  are  interested  in  the  case  where  <b.  is  the  incident  field 

xnc 

excited  by  the  source  in  the  exterior  region.  That  is, 
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(B19) 


9 


2  2 
(V  -  Y  > 


|>E 

me 


=  f. 


We  are  interested  in  the  case  where  4>ID  are  interior  volume 

n 

Dirichlet  eigenfunctions  satisfying  the  equation 


<n2  D  ,  .ID 
(V  -  YIn)  *n 


=  0 


(B20) 


and  the  Dirichlet  condition  on  the  surface 

4>*D  =  0  (B21) 

Substituting  Equation  (B19)  and  (B20)  into  Equation  (B18) , 
integrating  over  the  interior  volume,  applying  the  divergence 
theorem,  and  the  Dirichlet  boundary  condition,  yields 


/ 


.E  - 
incnI 


vof  dS 


-Kn-  V2)  / 


4>IDfE  dv 
n 


(B22 ) 


E 

Recalling  that  f  is  nonzero  only  in  V_  ,  we  rewrite 
Equation  (B22)  as 


/(-VO*Le«-  (l2- yll)  f  (B231 

S  Vj. 

In  order  to  draw  our  desired  conclusions  from  (B23)  without 
leaving  any  gaps,  we  would  have  to  reproduce  the  results 
connecting  the  surface  and  volume  solutions  for  Neumann 
boundary  conditions  contained  in  Section  IV  all  over  again 
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for  the  Dirichlet  boundary  conditions.  At  this  point  we  will 
not  provide  these  details  and  we  will  thus  have  the  indicated 
gap;  however,  we  see  no  difficulty  in  eliminating  this  gap. 
Assuming  the  Dirichlet  connection  yields  the  same  results  as 
the  Neumann  connection,  it  then  follows  that  the  volume  eigen¬ 
value  Yx  is  a  zero  of  the  surface  eigenvalue  XI,S'D. 

In  Dn 

According  to  Equation  (B15) ,  it  follows  that  y  is  an  extra- 

E  s 

neous  zero  of  A  '  which  we  can  formally  state  as 


SE,S,EX  .  cyD 


Im 


(B24 ) 


where  c  is  the  acoustic  propagation  velocity  as  opposed  to 

the  quantity  defined  by  Equation  (B2),  and  m  can  represent 

a  multiple  counting  index.  To  proceed  with  the  argument  it 

is  necessary  to  relate  -n^  •  V<j>^E  and  the  which  satisfy  (E 

A  Dirichlet  extension  of  the  Neumann  arguments  presented  in 

Section  IV  would  be  necessary  to  conclude  that  a*(yEn)  = 

A  •  V<tID  where  A  is  just  a  proportionality  constant, 
n  I  n  n 

Assuming  the  described  Dirichlet  equivalents  of  the  Neumann 
arguments  have  been  established,  it  follows  from  (B23)  that 


f  oI(sE'?'i:X)*=  (sE,!f'EX)dS  =  0  (B25) 

J  n '  nn’  c'  nn*  ' 

S 


This  can  be  written  as 


/  I*  E,s,EX.  E  ,E,s,EX.\ 
l°n  (snn*  ♦inc(snn'  ])  ~ 


and  using  Equation  (B16) ,  we  write  this  as 


I. Et  E,s,EX.  E  ,_E,S,EX  \  _ 

\  n  (snn'  } '  ^inc(snn'  )/  "  0 


(B26) 


(B27) 
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I 


The  left-hand  side  of  Equation  (B27)  is  the  original  defini- 
£  E  s  EX 

tion  N  ( s  )  and  we  now  have  our  desired  result, 

s  nn 

Having  obtained  the  desired  result,  we  could  end  this 
appendix.  Instead,  we  will  present  a  number  of  results  that 
have  some  aspects  in  common  with  the  results  that  have  been 
presented.  To  do  this,  we  summarize  the  essential  forms  of 
the  scalar  operators  that  have  already  been  defined.  For 
the  external  Neumann  problem,  we  have 

Lg  »  |  -  K  ,  (B28) 

for  the  internal  Neumann  problem  we  have 

Lg  =  |  +  K  ,  (B29) 

for  the  internal  Dirichlet  problem  we  have 

lID  “  \  ~  K+*  ,  (B30) 

and  for  the  external  Dirichlet  problem,  it  can  be  shown  that 


L__  = 

1  t* 

=■  +  K 

(B31 ) 

ED 

2 

From  these  forms  the  set  equalities  follow 

K} 

-  K> 

(B32) 

{  °n  } 

-  H> 

(B33) 
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K+}  -  K+}  <B34> 

K+}  ”  K+}  IB35> 

K’VKi 

k+*}  -  ki  ib37> 

{♦r*}  -  {°n }  (B3SI 

{Xn'S}  -  K'S'°}  (B39> 

as  well  as  Equations  (B15)  and  (B16) ,  which  played  a  central 
role  in  this  appendix.  In  addition,  the  following  eigenvalue 
equalities  hold. 


Xn'S  +  Xn'S  *  1  <B4(» 

Xn'S'D  +  ^'S'D  =  1  (B41) 

n  n 

Xn'S  +  =  1  (B42) 

n  n 

Xi'S  +  ll'S,D  =  1  (B43) 

n  n 

Finally,  for  the  spheroid,  which  includes  the  sphere  as  a 
special  case,  we  have 

K+*  =  K  (B44) 
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and  we  have  the  following  additional  set  equalities: 


K}  - 

un 

(B45) 

K>  - 

{♦n 

(B46) 

laEtn 

(B47 ) 

l  n  f 

\  n  f 

(aM  = 

(B48) 

l  n  f 

l  n  f 

For  the  sphere,  the  eigenfunctions  can  be  chosen  to  be  real, 
thus  eliminating  the  need  for  the  conjugate.  This  aspect 
of  the  sphere  solution  is  sufficient  to  show  that  K  and  K+ 
commute,  thus  making  K  normal. 
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APPENDIX  C.  A  USEFUL  SCALAR  SEM  RELATIONSHIP 


In  this  appendix  we  will  present  an  alternate  represen- 
ct  s 

tation  for  an^,  which  plays  a  prominent  role  in  the  final  SEM 
expansion  given  by  Equation  (60) .  There  are  several  alter¬ 
nate  forms  already  possible  for  this  quantity,  based  on 
material  already  presented  in  this  report,  and  the  form  we 
now  choose  is 


a 


a ,  s 
nn' 


(♦?«•> 


Tmc 


(s)) 


a,s 


n 


s=s 


nn 


s 

I 


(Cl) 


Implied  by  this  form  is  a  normalization  of  the  adjoint  eigen¬ 
function  that  appears  in  the  expression  as  well  as  a  normal¬ 
ization  of  eigenfunction  <f>^.  The  normalization  that  has  been 
implied  so  far  in  the  text  is 

(*>),  4>“(s>)  =  1  (C2) 


Such  a  normalization  was  always  possible  because  we  could 

have  chosen  the  unnormalized  functions  <J>a^  and  and  formed 
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The  normalized  eigenfunction  and  its  adjoint  can  then  be 
obtained  by  dividing  each  by  (Q®)^^.  We  are  now  in  a  posi¬ 
tion  to  express  Equation  (60)  in  terms  of  the  unnormalized 
eigenfunctions  as 
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where 
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It  is  the  denominator  of  the  last  expression  that  can  be 
written  in  an  alternate  form.  To  accomplish  this,  we  write 
the  original  eigenvalue  equation 


L“(s)^Q(s,r) 
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and  take  the  derivative  of  both  sides  with  respect  to  s  to 
obtain 
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Next  we  take  the  inner  product  of  all  terms  with  (f)^  to  obtain 


(C 


3Ls 

3s 


nQ 


.•!Sau 

Js  3s  } 


dX 


a ,  s 


ds 


a 

Q  +  X 
n 


a ,  s  /  a ,  + 
n  rnQ 


3<j> 


nQ' 


3s 


(C8) 


80 


The  second  term  can  be  written  as 


ha+ 

nQ 


Substituting  Equation  (C9)  into  Equation  (C8)  and  canceling 
terms  yields 
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The  desired  expression  is  obtained  by  substituting  (CIO)  into 
(C6)  and  it  is 
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We  now  have  a  choice  of  dealing  with  an  explicit  normalization 
procedure  implied  by  Equation  (C6)  or  the  normalization  insen¬ 
sitive  procedure  implied  by  Equation  (Cll) .  It  should  be 
noted  that  the  electromagnetic  SEM  equivalent  of  the  material 
presented  in  this  appendix  .has  long  been  established  using 
virtually  the  same  procedure. 
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APPENDIX  D.  A  TIME  FOR  THE  EQUIVALENCE  OF  CLASS  1 
AND  CLASS  2  COUPLING  COEFFICIENTS 

In  this  appendix  we  show  that  when  a  perfectly  conducting 
body  is  illuminated  by  an  electromagnetic  pulse  of  finite 
width  T  the  SEM  responses  with  class  1  and  class  2  coupling 
coefficients  are  identical  for  times  t  ^  T  +  (L/c)  +  t  , 
where  L  is  the  maximum  body  dimension  in  the  direction  of  prop¬ 
agation  and  t=tQ  corresponds  to  the  instant  at  which  the  pulse 
wavefront  first  hits  the  body. 

First  we  show  that  this  is  true  for  a  delta  function  inci¬ 
dent  pulse,  i.e.,  for  T=0.  The  coupling  coefficients  are 
defined  as  follows 

=  Ae(Ya“Y)ctoy*jlnc  (rfya)  .  $a(r)dS  (Dl) 
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where  A  is  a  common  factor  and  $  (r)  is  a  coupling  vector. 
For  a  delta  function  pulse 
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where,  for  simplicity,  the  direction  of  propagation  has  been 
chosen  along  the  z  axis,  with  z=0  at  the  "center"  of  the  body, 
p  =  n  x  h,  n  is  the  unit  normal  to  the  surface  of  the  body  and 
h  is  the  polarization  vector  for  the  incident  magnetic  field. 
In  view  of  Equation  (D3) ,  Equation  (D2)  can  be  rewritten  as 
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where 
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If  we  Laplace  invert  the  SEM  simple-pole  expansion 
term  is 
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where  C„  is  the  Bromwich  path.  For  times  t> (z/c) 
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The  pulse  wavefront  has  just  passed  the  body  at  t=(L/c)  +  tQ 
and 
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The  corresponding  term  for  class  1  is 
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and  is  identical  to  T*2'  in  Equation  (D5) 

For  a  general  incident  pulse 
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We  see  that  the  above  expressions  involve  the  product  of  £(y) 
with  another  function  of  y  and,  consequently,  the  convolution 
theorem  can  be  employed.  If  we  rewrite  Tv  '  as 
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then 
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where  F(t)  and  P(t)  are  the  inverse  Laplace  transforms  of  f(y) 
and  p (y)  respectively.  Similarly,  for  T*1^  we  can  write 
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The  calculation  for  the  delta  function  pulse  showed  that  P(t) 
is  equal  to  Q(t)  for  t  >  (L/c)  +  t  .  We  will  now  show  that 

(2)  .  °  /i  \ 

T'  given  by  Equation  (  6)  is  equal  to  Tv/  given  by  Equation 
(D7)  for  t  (L/c)  +  T  +  tQ  where  T  is  the  pulse  width  (Fig.  Dl)  . 
This  is  done  simply  by  plotting  the  products  F(t-x)  P(x), 

F(t-x)  Q(x),  and  observing  with  the  aid  of  Figure  Dl  that 
for  t  >  (L/c)  +  T  +  t  ,  and  T ^  are  indeed  equal. 

O 

Notice  that  as  T-*»  the  SEM  responses  with  class  1  and  class  2 
coupling  coefficients  are  different  for  all  times  as  noted 
in  Reference  1. 
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/  t- T  T  t 

(L/c)  +  tQ 


Figure  Dl.  The  four  graphs  show  how  to  plot  the  products 
F(t-x)  P(t)  and  F (t-T )  Q(t). 
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APPENDIX  E.  MAGNETOSTATIC  IMPLICATIONS  OF 

COUPLING  COEFFICIENTS 

In  this  appendix  we  calculate  the  static  current  density 
induced  on  the  surface  of  a  perfectly  conducting  sphere  immersed 
in  a  homogeneous  magnetic  field  by  calculating  the  class  2 
SEM  current  density  for  a  step  function  incident  pulse  as  t-*«>. 
This  static  response  has  the  correct  6  and  <j>  dependence,  as 
it  should,  but  an  incorrect  numerical  coefficient. 

To  facilitate  the  analysis  we  assume  that  the  direction 

inc 

of  propagation  is  in  the  z  direction  and  that  H  is  in  the 
y  direction.  This  situation  corresponds  to  a  polarization 
index  p=2  with  e^e^,  e2=ex,  ez=ey  {ex=0,  4>1=n)  (Ref.  12).  Thus, 
following  Baum's  notation  and  the  results  in  Reference  12 


where  o,e  stand  for  odd,  even  and  only  the  m=l  terms  are  non¬ 
zero.  The  class  2  coupling  coefficients  have  been  evaluated 
in  Reference  1. 
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where  f(s)  is  the  Laplace  transform  of  f(t).  Thus  we  should 
evaluate  cyJ(r,Y)  as  y+o.  We  have 
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We  have 


(a/c>  Dl,l,o  =  1  ^(Qa  -  -1 

i,(-l)  =  -1,(1)= -.36788,  i2(-l)  =  i2(l)  =  .07156 

R,  =  =  -cos  4>  +cos  0  sin  <p 

—1 ,1,0  o  <p 
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and 


lim  J(r,t)  =  .45568(-  cos  <J>  eQ  +  cos  0  sin  4>  §.) 
t-*°°  °  *** 

The  correct  response  has  a  numerical  coefficient  equal  to 
1.5  instead  of  0.45568.  Since  the  SEM  solution  which 
employs  class  1  coupling  coefficients  is  equivalent  to  the 
Mie  solution,  the  use  of  class  1  coupling  coefficients  must 
yield  the  correct  magnetostatic  solution. 
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